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developed and its use was demonstrated through a design 
example. "SPRMD" was written in FORTRAN for use on an IBM 
PC/AT. It combined several existing codes ("MICROPEP", 
"GRAINS", "ROCKET", etc.) and used the performance loss 
estimation methods suggested by the AGARD Propulsion and 
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IT. INTERODUCT IGN 


A. BACKGROUND 

Solid propellant rocket motor design is a very complex 
process which often requires many iterations in order to 
accomplish a successful design. Usually, the first step in 
the design process it to obtain accurate input data oon mis- 
Sile system restraints. Examples of these restraintS are 
ballistic requirements (which can include thrust, burning 
time, total impulse) and mission/vehicle requirements (which 
can include motor length, diameter, allowed hazard and plume 
properties, altitude, etc.). These parameters are often re- 
ferred to as independent parameters. Through the use of 
these independent parameters, the required propellant burn 
time and properties can be roughly estimated. Then, theo- 
retical motor performance and losses can be estimated. The 
thrust-time requirements together with the expected struc- 
tural loads permit initial selection of a grain design. 

The next step of solid propellant rocket motor design is 
to use the propellant properties and grain design to preuict 
the pressure- and thrust-time profiles a the total impulse 
of the motor. Propellant properties and grain design are 
then iterated until these independent variable requirements 
are met, together with acceptable propellant stress levels. 


A final step in the missile design process is to use the 
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thrust-time profile in an external aerodynamics code _ to 
ensure that range and terminal velocity requirements are 
met. 

During the past decade, there has been significant pro- 
gress in the automation of solid rocket motor design. It is 
mainly because of the extensive use of digital computers, 
which greatly facilitate the iterative type computations 
mentioned above (Ref. 1). In recent years, many of the 
design processes have been accomplished using separate com- 
puter codes, 1.@e., adiabatic equilibrium combustion codes, 
codes which estimate performance losses, grain geometry 
optimization codes, etc. These codes have, in some cases, 
been combined into much larger design codes, and they are 
often difficult and expensive to run. Recently, PC versions 
of some rocket motor design codes have become available. 
However, they must still be used individually, and were not 
all available for use on the same computer. The use of 
these codes was also very inconvenient in that the same 
parameters often had different variable names in the various 
codes. The user had to transfer or convert the data from 
one code to the next. 

In addition, a recent AGARD publication has provided 
recommended emirical equations which can be used _ for rapid 
and accurate estimation of pverformance losses for aluminized 


propellants. 
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Bs OBGEeT IVES 

The first objective was to outline a suitable "Solid 
Propellant Rocket Motor Preliminary Design Procedure" which, 
of course, 1S not a unique procedure. The second objective 
was to create, modify, and/or combine the PC versions of the 
necessary design codes into one interactive code which could 
be conveniently used for preliminary design on an IBM PC/AT 
microcomputer. The final objective was to demonstrate the 
program by designing a tactical solid propellant rocket 


motor combustor. 


C. APPROACH 

The approach to achieve the objectives was first to 
detail an acceptable preliminary design process. The algo- 
rithm of the design process was then followed, utilizing 
several existing codes and loss estimation methods. as 
appropriate, to obtain an integrated design code for use on 
the IBM PC/AT. The programs to be used were (1) “MICROPEP" 
(a microcomputer version of a pronellant evaluation program 
from China Lake, Naval Weapons Center), which calculates the 
adiabatic, equilibrium combustion process, and theoretical 
motor performance; (2) "GRAINS" (a radial spoked star grain 
design program from CSD); (3) "ROCKET (a Lockheed internal 
ballistics program), which is a one-dimensional flow code 
which predicts pressure and thrust versuS time; and (4) 


"FLYIT” (a flight simulation program from China Lake, Naval 
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Weapons Center), which performs simplified trajectory calcu- 
lations. In addition, a code based upon AGARD-AR-230 was 
written for estimation of performance losses for aluminized 
propellants. Capabilities of the final program were then to 
be demonstrated by considering the design of a tactical 
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Ii. SOLID PROPELLANT ROCKED MOTOR DESTGNSEROGESS 


A. GENERAL 

The method of solid rocket motor design adopted in this 
work is shown in Figure 1. It was assumed that the required 
information on fundamental propellant properties (burning 
rates, temperature sensitivity, etc.) was known. In 
addition, the design process discussed herein did not 
include any stress analysis for the propellant grain or 


Case. 


B. INPUTS 
The inputs that are required to satisfy mission objec- 
tives (the independent parameters) were as follows: 
1. Ballistics Performance Requirements 
a. Average thrust (F) or thrust-time profiles, 
Bb. Burn) times Geos 
c. Total impulse (I;)- 


These parameters are interrelated as shown in 


eb 
Figure 2. Thus I; = f Fdt (2) 
O 


2. Mission/Vehicle Constraints 


These requirements are directly related to the motor 


design and include: 
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DESIGN PROCESS 


i 


a. Envelope 
Tne allowable envelope (motor lenath L, motor 
inside diameter D) is the fundamental constraint on the 
motor geometry, which can be |. used to estimate the required 
burning rate (r) and total volume available (V,) for the 


propellant by following relations: 

r = (D/4)/tp ( 28e5 
(1.@., a web fraction of approximately 5) 

Va = (D*/4) L (2.45 


b. Maximum Expected Operating Pressure (MEOP) 

The MEOP is set by the structural limit of the 
motor case. The maximum operating pressure at maximum pro- 
pellant temperature (Pmaxtp) is less than the MEOP by the 
factor of reproducibility tolerance (see Figure 2). Pmaxtp 


can be initially estimated from the following equation 
Pmaxtp=MEOP (l-reproducibility tolerance) (2 .5® 


c. Temperature Range and Altitude 
The operating propellant temperature range 
effects the selection of the propellant through the required 
burning rate temperature sensitivity. The design altitude, 
together with the nominal chamber pressure determine the 


exhaust nozzle area ratio. 
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FIGURE 2. 


TYPICAL PRESSURE OR THRUST VERSUS TIME PROFILES 


Le 


d. Allowed Hazard Rating and Allowed Exhaust 
Signature 


The hazard rating requirements have to do with 
whether the propellant 1s detonable or non-detonable. 
Modern high energy propellants are often detonable. The 
allowed exhaust signature affects the propellant selection 


through requirements on exhaust smoke and radiation. 


C. PROPELLANT SELECTION 

From the ballistics performance and mission/vehicile con- 
straints, it is possible to make a preliminary selection of 
the propellant type that is most likely to give the required 
performance, internal ballistics, flame temperature and 
mechanical properties, as well as the necessary storage 
Stability, the best hazard properties and exhaust signature 
(REE. 92). A particular propellant will have known (or 
estimated) burning rate (r), temperature sensitivity (tx), 
density (Pp), and composition. Then the nominal chamber 
Dr eseiiee Ge. ) can be calculated by usingethe» sequaieien 


nom 
(Ref. 3): 


Po = Pmaxtp @XD (™K (Tnom - Tmax)? (2.6) 
nom 


D. THEORETICAL PERFORMANCE COMPUTATION 
Following the propellant selection, theoretical per- 
formance computations can be made for the specified chamber 


and ambient pressures. A program such as PEPCODE (Ref. 4), 
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or the PC version "MICROPEP", will provide data such as 
thrust coefficient, characteristic exhaust velocity, shift- 
ing and frozen equilibrium performance results, optimum noz- 
zle expansion ratio, mole fraction of condensed species in 


the exhaust, etc. 


E. NOZZLE GEOMETRY AND MATERIAL SELECTION 

The information on theoretical combustion temperature, 
pressure, gas composition, and burn time, together with 
geometry constraints, permits the nozzle configuration and 


meeemeial to be selected. 


Fe. PERFORMANCE LOSS ESTIMATION 
It was assumed that the delivered specific impulse could 


be written: 


iep = 7) 


SD NC Nc 
where aD 1s the theoretical specific impulse calculated 
by “MICROPEP"; and ore and nc* are the thrust coefficient 
efficiency and characteristic velocity efficiency (or com- 
bustion efficiency), respectively. 

NC - was computed by summing the effects from the 
following losses 

1. Divergence losses (epry),s 


2. Two-phase flow losses (ep), 


3. Boundary layer loss (epy), 


he 


4. Kaneticess losee( cami 
5. Submergence loss (esyp), 
6. Erosion loss (eprpos) - 


Thus, NC ¢ = 1l-(epnzyteppt epy.t ex int Esupt EpRog)/100 (2 Joe 


The combusion efficiency depends primarily on the 
residence time (t,;), and the latter can be expressed as 
ty = (Vol/m) (P, M/12 R Tg) ( 2h 
where 
- Vol i che instantaneous volume of the combusor (in), 
= wee is the chamber pressure (lbf/in?2), 
- T, is the chamber temperature (9R), 
- MM is the molecular wero, 


- R is the universal gas constant (ft lbf/lb,,51,°R),and 


| 
ge 


is the mass discharge rate (lbm/sec). 

For example, Figure 3 (Ref. 2) is a plot of specific impulse 
versuS average residence time for aluminized propellants. 
It shows that the residence time must be somewhere between 
10 to 15 msec to obtain a combustion efficiency between 90% 
and 95%. Increasing the free volume of the combustion 
chamber will improve the residence time. In this work the 
initial volume was used in the residence time calculation, 


1.e., the worst case. 


G. THROAT AREA AND EXIT AREA 
Accurately sizing the throat and exit areas is one of 


the major goals of the solid propellant rocket motor design. 
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FIGURE 3. 


EFFECT OF RESIDENCE TIME ON Isp EFFICIENCY (Ref 2) 
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These areas are dependent on several other variables, 
such as thrust coefficient efficiency, combustion 
efficiency, specific impulse, and etc.. In the sizing 
process, it is convenient to start by using the theoretical 
values. The effects of losses can then be handled itera- 
tively until the solution converges to the final areas. The 
rest of this section will present the basic algorithm for 
the throat and exit area sizing process. 

The basic definition of specific impulse can _ be ex- 


pressed as (Ref. 3): 


ee = F/m (2 
where 

F = Ce, PcAth (2.1) 
Then 

Ath = E (2 «lee 


Cfth "Cg Pe 


The theoretical throat area can be estimated by using 


h 
theoretical performance computation). The losses which are 


fo for shifting equilibrium flow (from "MICROPEP" — ae 


calculated to wovotain Eee depend upon the output from 
"MICROPEP" and the throat area. Thus, iteration is required 
to obtain both the losses and throat area. 

In summary, the calculation can be conducted in the 


Following steps: 


Ae 


Pccluetmearems weusing equation” (2.11). 
2. Calculate nc_ as discussed above. 
fe 
3. Calculate Arp Wsmngeeatuatlon (2.12). 
4. Repeat steps 2. and 3. until the values for throat 
area and ier converge. 
In the current program the iteration was terminated when 
there was less than a 1% change in Ayy, Per iteration. 
The exit area can then be calculated from 
Ae = Ath & (Zero) 
where € is the exit to throat area ratio. e¢ is calculated 
in "MICROPEP",, assuming that the nozzle exit pressure 1S 
equal to the local ambient pressure (one-dimensional, ideal 


flow is assumed). 


H. GRAIN DESIGN 

The primary purpose of grain design is to shape the pro- 
pellant grain such that the rocket motor can deliver the 
required performance (thrust-time or pressure-time profile). 

In order to be able to select the general configuration 
of the grain geometry, several parameters must first be 
determined: the web fraction (W¢), volumetric loading frac- 
tion Gee andmernc length>to-diameter ratiog(L/D) ~(Ref. 5). 
L/D 1s known from the input data. Web fraction is the frac- 
tion of the motor radius that is filled with propellant. 
Thus, 


rth 


(b/2) (2 eid) 


Wr = 


ZS 


The propellant burning rate (c) is generally expressed in 


the form: 


n 

r= apc (2. ay 
where a and n are empirical constants. They are assumed to 
be known for the selected propellant. The volumetric load- 


ing fraction is defined as the ratio of the volume occupied 
by the propellant (Vp) to the total volume available in the 


motor for propellant (Va). Thus, 








Vp 
where 
nD 
Va= (=,-) (2 aaa 


The required propellant volume is found as shown below 


a W 
le = Step Lestes = eee - th =e leans (2. 1s 


Thus, Wp = It/Isp (2 le 


This yields the required propellant weight (Wy) and 


Vp = (2.20) 


Shs 


where Pn is the density or specific weight of the selected 
propellant. With We, V2, and L/D known, an initial selec- 
tion of grain geometry can be made (for example, see Ref. 
Sie For a selected grain configuration, the web versus 


burning surface area (which is the data needed for 
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thrust and pressure versus time calculations) Caneaoe 


computed. 


bh. THRUST AND PRESSURE VERSUS TIME CALCULATIONS 

The thrust and pressure versus time calculations are 
made using conservation of mass for the combustor(Ref. 3). 
The mass addition from the burning propellant (Mg) 1s 
equal to the sum of the mass accumulated in the combustor 
Migeene Mass exhausted from the nozzle (my) as shown in 


Figure 4. Thus, 


Mg = dM/dt + mp (22a) 


Mg is the mass of gas generated from the burning propellant. 
® > Nn 
Mg = )PpApr = )ppApaPc (C22) 
i 


The summation is over all propellant grains that are burning 


at the same time. 

















= 1s the mass accumulating in the combustor. 
= —<—( pgVol) = og te oul cfg 
where : (2) 
. d 1 dP 
=. LrAp ae “a * REG dE 
then 
RT¢ RT 6 
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FIGURE 4. | 


MASS CONSERVATION IN A SOLID PROPELLANT MOTOR 
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Mn is the mass exhausted from the nozzle. For chocked flow, 








DP. A 
as cae oc 
fin = SERS (226) 
Combining equations 2.21 through 2.25 results in 
P ap PeAL eG 

FopApaP2 = = (—-S) Wohi Yotj se + (Ss (DQE2S) 

1 RTo e 
selving for dP/dt 

dPo 1 \r n _ PcAthIc dVol 

=n (Fey (RT) PpApDpaPG - —er—) Pol-ae ) | (Perez) ) 


Peet iene 2627 can be numerically integrated to give P,(t). 
The thrust-time behavior can then be obtained from 


F(t) = ng. Cg(t)Po(t)Agn(t) (228) 


J. TERMINAL VELOCITY AND RANGE CALCULATIONS 

The last step of solid propellant rocket motor design is 
to incorporate the designed motor into the desired airframe 
and make a flight simulation. The purpose of this flight 
Simulation is to verify that the motor can produce the 


desired range and terminal velocity. 
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III.  MECRGEEE 


"MICROPEP" is a microcomputer version of the propellant 
Evaluation Program (PEP), published in Naval Weapons Center 
Report NWC-TP-6037, "Theoretical Computations of Equilibrium 
Compositions, Thermodynamic Properties, and Performance 


Characteristics of Propellant Systems," by D. R. Cruise, 
April 1979 (Ref. 4). The code (MICROPEP) can handle 12 dif- 
ferent input ingredients and 200 combustion species. Te 
performs the calculations for the thermodynamic properties 
and performance of propellant systems. 
A. BASIC ASSUMPTIONS 

The basic assumptions of "MICROPEP" are as follows: 

1. JE QUaaIE i. tum adiabatic combustion at the nozzle 
entrance. 

2. One-dimensional, isentropic expansion in the nozzle in 
which the flow can be either a shifting equilibrium 
flow (assuming chemical equilibrium 1s _ maintained 
throughout the expansion process) or frozen composi- 
tion flow (assuming the chemical composition of the 
the flow remains the same throughout the expansion 
process). 

3. Gases were considered to be perfect. 

4. Nozzle exit pressure is assumed to be equal to the 


ambient pressure. 
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B. INPUTS 
Inputs required by "MICROPEP' are: 

1. Composition and mass of each ingredient. Normally the 
weights of ingredients are chosen to add up to 100 gq. 
NOTE: Ingredients may be input by the user or they 

may be selected from an existing data file. In 
the latter case ingredient serial numbers are 
used. 

For assistance, use README, which should be on 
tne PEPCODED.DAF Data File. 

2. Combustion pressure. 

Be ee. pressure. This defines the limits for the 


expansion process. 


Geeoo re UTS 
The program calculates: 

me Combustion product COnposTrlon, properties , and 
temperature at the nozzle entrance. 

m™ wNozzre exit PleGduick  COMmpesition, properties, and 
temperature for two cases: 
a. Steen equilibrium flow 
b. Frozen composition flow. 

3. Required throat area and exit area with no losses, for 
the same two cases. 

4. Theoretical specific impulse, thrust coefficient, and 
characteristic exhaust velocity for both shifting 
equilibrium and frozen composition flow. 


ao 





IV. LOSS MECHANISMS 


For estimation of the performance losses, a set of 
empirical formulas Suggested by the Propulsion and 
Energetics Panel of AGARD in AGARD-AR-230 (Ref. 6) were 
used. Similar methods can also be found in NASA Report SP- 
8039 (Ref. 7) "Solid Rocket Motor Performance Analysis and 
Prediction." It is assumed, as in equation 2.7, that the 
combustion efficiency (nc*) and thrust coefficient effici- 
SN yes )) Selo iS Wiwesciceel is independently. The empirical 


formulas are valid only for aluminized propellants. 


A. THRUST CORBETGGENT ERR TRGWEri ete 
The thrust coefficient efficiency ce? can be expressed 


as in equation 2.8 


Ne, ~ TO VEDIV! TP SBME KIN 9S UB \ieme neue me 


Empirical equations for each of the losses (e's) will be 
presented and briefly discussed below. More details can be 
found in Reference 6. 

1. Divergence Loss (DIV) 

This 1s the loss due to the radial velocity com- 
ponent of the gas at the nozzle exit plane 


EpiIy = 90 [1-Cos Sanne) (4.1) 
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a is the nozzle half angle and dex 1s the exit angle for a 
contoured nozzle. For a contoured nozzle, a 1s measured to 
a line drawn from the exit and tangent to the wall at the 
throat. 

2. Kinetics Loss (*KIN) 

This is the loss (or reduction in performance from 
the value of shifting equilibrium flow) associated with not 
attaining chemical equilibrium throughout the nozzle expan- 
S1lon process. It 1S estimated to be 1/3 of the fractional 
difference between the specific impulse for shifting equi- 


librium (Igp ) and frozen composition (Isp .)- Thus, 
S, 


I 
oon (|) = © (4.2) 


hn I 
SPs 


Both Wee and Isp-¢ are ebtainead Erom ‘MICROPEP.” 


3. Boundary Layer Loss (*BL) 


The boundary layer loss (ep,) 1S expressed as 


P 0.8 Pp e 

Cc &: 
EBL = Cy ——prz [1+2 exp(-Cz G5 

Dth Dth 


ty) ] [1+0.016(e-9) ] 


(4.3) 


This loss accounts for both reduced flow areas and transient 
heat loss. In equation 4.3, P, is in psi, Dy, is in inches. 


The C, and Co coefficients are 


Sl 


for ordinary nozzles Cy aqme 0 «00265 
Co = 0.000937 
for steel nozzles Gy = O00 0506 
C> = Une 
4. Two-Phase Flow Loss (TP) 

This loss is due primarily to the velocity lag of 
the condensed species of the combustion products as_ the 
mixture passes through the exhaust nozzle, 

epp = C3 a Dp's (4.4) 

pet 0.08 Depac 
where 

a. De, Us the diametereot theameineinenee 

be Dp is the mean A150, part ieren diame ten in 

microns and is calculated from 
Dp = 2 iG Dp p02 4692 
ec. PP. 1S chamber prescimees ingore 
d. C's are dependent on the diameter of the throat 
Cn = 0.5 
Dyes ee oe Catia So Ce = 
it < Dth < Zs See eee Cs = l Ce = Ona 
Dep? 2 and D,<4: C3 = 13.4 Cy = 0.8 Ce = Oe 
Depo 
4<Dp <8: C3 = 10.2 Ce = 0.8 Ce — ae 
D+ p72 and D,>8: C3 = 7.98 Ce = 0.8 Ce = (Ge 
e. Mpg is the mole fraction of consdensed phase in 


moles/100 gm of 


3 


Z 


reactants. 





5. Submergence Loss (*SUB) 


This loss is expressed as 


age ) me? @4.5) 


ESUB = 0.0684 ( 





It is due to the nozzle configuration for a submerged nozzle 
where S is the submergence length in inches (the length that 
the nozzle is imbedded in the combustion chamber). 
6. Nozzle Erosion Loss (*EROS) 
This loss results from nozzle erosion during the 


burn and is expressed as 


, i 


Tsp. y, is the theoretical Isp for the initial nozzle expan- 
Sion ratio and toe is the theoretical Isp for the mean 
expansion ratio. I spn is determined by interpolation of 


Soe Vemsus. ele nozzle expansion ratio data from 


SUeROPEP. ” 


Be. COMBUSTION EFFICIENCY 

NASA Report SP-8064 (Ref. 2) states that the combustion 
Seercleieine+) Is .d@termined by “the completeness of metal 
combustion within the motor and by the degree to which com- 
bustion products reach chemical equilibrium among 
themselves." . These effects depend primarily on the 


residence time (tr) as expressed in equation 2.9. An nex of 
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90-95% can be attained if residence time is 10-15 msec (see 


Figure 3). 


C. PROGRAM 
The subroutines which solve for the thrust and combus- 
tion efficiencies are listed in Appendix A, Subroutine PFE 
and Subroutine LOSCF. Basically, they followed the method 
discussed in Chapter II. The algorithm is shown below 
1. Initialize Cre (coefficient of thrust) as Cf, 
(theoretical coefficient of thrust from "MICROPEP" 
FOr Shishelne, COU: Lor lume 


2- Calculate Arp and Dep? 


7 ee 


7 Ee 
Dth Ey (cE Bot 


3. Calculate ce using the relation 


ies I-( epnry+epyt eppt ex nt esupt epros) /100 


4. Caiculate Cr using 


Cea=ec n 
f Eih Ce 
5 . | VCOEECCES Dar mu GrnG 
= eZ 
Deh = 2 (EL) 
Ce ° Pat 
6. Repeat Steps 3 to 5 until Ay, converges. 
The reason that Ayp, was used for the iteration variable was 


that no (through ep,, epp, and Ecijin) depemds upon Ath: 
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V. _ GRAINS 


"GRAINS" iS PC based code for the design of radial sided 
Star grains. The original version waS written by G. J. 
Woten and R. McCormick of CSD, United Technologies, November 
1984. It was then modified by LT G. Liston of AFWAL/Port, 
December 1984. The general shape of this grain and the 
geometric variables are shown in Figure 5. 

"GRAINS" can be used in either of two modes: the 
"Design Mode" or the "Burnback Model." 

A. “DESIGN MODE" 

In the "Design mode" the program solves for geometries 
Petemsaticoty input requirements on the number of star points 
required port area, perimeter factor, web, ane: CbC%; 
together with the allowed tolerances on these parameters. 

The geometric constraints for the grain are determined 
from the propellant characteristic and motor constraints. 
For example: 


Ft b 


)| 


Aport = (7)p4[(1 - Sar 


D, tp and F are motor inputs requirements and I,, and Pp are 


D 


determined once the propellant has been selected. 


3D 





FIGURE 5. 


NOMENCLATURE FOR "GRAINS" DESIGN MODE 
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B. "BURNBACK MODE" 

The "burnback mode" uses the length of the _ grain, 
initial port area, perimeter, and the results from "Design 
Mode" to calculate the perimeter, perimeter factor, port 


area, and burning surface area as functions of the web 


burned. 


C. COMMENTS ON PROGRAM OPERATION 

Care must be used for the inputs and input tolerances, 
otherwise no acceptable solutions will be obtained. It is 
best to start with large allowable tolerances on Aport, etc; 
mevoncally 20%). ieee Weel ec Oluit Ons are. foundeain 
"Design" for the number of star points, range of web 
thicknesses, etc., the range of the input variables can be 


reduced along with the allowable tolerances. 


oe 


VE. ROCKED 


"ROCKET" is also a PC-based computer code. It was 
Originally written by The Lockheed Company for a mainframe 
computer. Later, it was transferred for use on the IBM-PC 
by J. P. Francis of the Naval Weapons Center (August 1985). 
The function of "ROCKET" is to calculate the pressure-time 
and thrust-time profiles for the motor, including theoreti- 
cal calculations for both ignition and talloff transienees 
The description of the motor can include igniters, multiple 
propellant grains, and motor insulation. In this thesis 
work the program was restricted to amotor with a single 
grain configuration. 

A. BASIC ASSUMPTIONS 
The assumptions made in "ROCKET" are: 


1. The combustion products are ideal gases. 


2. Propellant burning rate follows the expression 
: n 
r= aPo - 


3. There are no effects of mass addition or erosive 
burning in the combustion chamber. 
A. C* Varies wit Po in the fommeot 


x Po 
CAE 
ref (p— 





ine 
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5. The temperature sensitivity (mg) 1S used to correct 
for the effects of propellant temperature on the 


propellant burning rate and chamber pressure. 


@ @ oc oO 
rp = L(70°) e0.0lop,(T, Qe) 
Dp 
where 
OH = l= ati 


wt 


B. ANALYTICAL BACKGROUND > 

Mires program 1s based upon the method discussed in 
Chapter II, 1.e., it numerically integrates equation (2.27) 
memeeeain Po(t). Them it calculates Ce with a fixed ratio 
of specific a (ey. Thrust is then calculated using 
equation 

F = Celt) Po (t)Agp(t) ACh 

where iA is nozzle divergence loss and Cp is the nozzle 
discharge coefficient. h and.Cp are only two parts of nem: 
In order to use ROCKET in the combined program, i was” set 


to unity and Cp was equated to ae 


ee iNePUTS 
Pietes Con ROGKET™ inclinde: 
l. General motor parameters: 
a. Throat area (A¢p) 
b. Nozzle exit area (Ag) 
c. Total motor volume (V,) 


d. Ambient temperature and pressure 


oo 


e. Nozzle throat opblug closure blowout pressure 
(nozzle is plugged until this pressure 1S 
exceeded, typically 35 psia) 

f. Initial pressure in the motor (Pzero, typically 
15 psia) 

g. Nozzle divergence loss A (set equal to 1 in this 
work ) 

he. Ratio of specific heats for chamber gases (y) 
(this was obtained from "MICROPEP" chamber 
conditions) 

i. Design pressure (P,) 

je Throat radial erosion rate (erosr). 

2. Motor element descriptions: 

a. Propellant burning rate (xr) 

b. Burning rate exponent (n) 

c. Propellant temperature sensitive (t,x) 

d. Burning rate reference pressure 

e. Propellant C* 

f. Pressure correction exponent for C* 

g.- Propellant density (Pp) 


h. A table of web versus hurning surface area. 


DD.) CUE is 
The outputs include: 
1. A pressure versus time table 


2. A thrust versus time table 
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3. The free volume, propellant consumed (weight), and 


throat area versus time table. 
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Vil. 


"FLYIT" is a three-degree-of-freedom (range, altitude 
and angle of attack) trajectory simulation program origi-~ 
nally developed by G. Burgner of the Propulsion Analysis 
Branch, Naval Weapons Center (NWC). A User's guide was 
later published by Y.G. Coenen (Ref. 8). The program can 
Simulate vertical-plane trajectories which include take-off, 
climbs, cruises, landing, and rocket boost. It 1s a tool to 
synthesize and analyze the trajectories of air-to-surface, 
and surface-to-surface missiles. 

The basic structure of "FLYIT” conSasts of five flies 
tional sub-groups, they are: 

- Atmospherics 

- Trajectory comtrois 

- Aerodynamics 

- Rocket boosters 

- Sustanier propulsion. 

Users are allowed to select a different option in each 
Functional. group to simulate different flight profiles with 
different propulsion systems. Through the use of "FLYIT" it 
is possible to estimate the terminal velocity and range of 
certain propulsion systems. 

"FLYIT" was not incorporated into the main body of the 


present program due to time limitations. 
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VIII. COMBINED PROGRAM 


The combined program ("SPRMD") is listed in Appendix 
and basically follows the design process discussed in 
Chapter II. The basic structure of this program 1s shown on 
Eeure 6. 

This program combines "MICROPEP" (to calculate the 
theoretical thermodynamic performance of the propellant), 
"PFE," and "LOSCF" (to calculate nc,, and to SSeimate oxy 
pemesmze the  ~throat and exit areas), "GRAINS" (to provide 
web burned versus burning surface area for radially sided 
star grains), and "ROCKET" (to compute the pressure-time 
and thrust-time OE ree). Inputs to the program include 
"INPUTO" (the initial mission requirement inputs), "INPUT1" 
Penouts for wi CcROBS> ), LNPUT2Z” (anputs for "GRAINS"), and 
meet iN’ (inputs for "ROCKET"). All of these programs com- 
prise an interactive program for solid propellant rocket 


motor preliminary design. 


A. ALGORITHM 
A simplified algorithm of this program can be briefly 
listed as follows: 
ie “INPUTO" 
Input the basic ballistic and mission/vehicle 


requirements. Calculate the required average 
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INPUT VARIABLES 







"MICROPEP" AT 2 P's oat. ge 


THEORETICAL PERFORMANCE 


AGARD LOSS CALCULATIONS 
= ESTIMATED ACTUAL PERFORMANCE 


“"GRAINS* 

Ab VS. WEB BURNED 

THRUST VS. TIME 
: PRESSURE US. TIME 

| | = 
[gj TERMINAL VELOCITY 

| ‘| RANGE 
FIGURE 6. 


FLOW CHART FOR THE COMBINED PROGRAM 
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requirements. Calculate the required average 

chamber pressure for a specific propellant = and 

propellant temperature. 

"“INPUT1" 

Input propellant composition, chamber pressure, and 

ambient pressure. Create a file for 'MICROPEP." 

“MEGROPEP™ 

Perform We hee performance calculations for 

the selected propellant. 

SPREE,” 

Performs the following functions: 

a. Estimates the combustion efficiency 

b. Calculates the delivered specific impulse 

ce. Sizes the throat area, exist area, calculates 
perimeter factor, perimeter, required web, 
required burning rate, and etc. 

EN PuUr2” 


Input the geometrical description of the star grain 


and propellant characteristics. Create a file to 
accommodate the processing of "GRAINS." Other grain 
design programs (such as _ SPP) can be readily 


incorporated. 
"GRAINS" 
Designs the star grain and calculates web versus 


burning surface area. 
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vm PRekK LEN 
Uses the results from the previous programs to create 
a file for "ROCKET." 

8 am ee ROCKER. 


Calculates the pressure-time and thrust-time profiles. 


Bz BRANCHING 

The algorithm discussed above also provides for the 
Capability of ~beanchenar For example, if the desired Po 
versus time 1s not attained, the user can either re-do the 


grain design and/or select other propellant properties. 
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IX. EXAMPLE OF DESIGN PROGRAM 


The initial requirements which are input into the pro- 
gram are shown in Table l, together with a selected value 
for the temperature sensitivity of the propellant (t,x). A 
reproducibility tolerance of 7% was chosen. 

The nominal operating pressure at Tp = 10°F Swasteaicu- 
lated to be 1030 psia and the required propellant burning 
rate for a web fraction of 0.5 was 0.615 in/sec. 

A composite propellant consisting of 78% AP, 12% HTPB, 
and 10% Al was selected. “MICROPEP" was then run for pres- 
sures of 250 psia and the nominal P,. value of 1030 psi. 


Outputs from "MICROPEP" are presented in Table 2. 
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TABLE 1 


INPUT PARAMETERS 


AVERAGE THRUST, LBF.. ccc cece ccc ccc 14250.0 
BURN TIME, SEC. . « os « « sisiemeletsie « « « « sreneneene 6.50 
MINIMUM THRUST REQUIRED, LBF.......--e. 8000.000000 
MOTOR LENGTH, INCH... cc cccccccvvccvcee 40.00 
MOTOR DIAMETER, INCH. cc we eee ccc ccc vcee 16.00 
MAXIMUM EXPECTED OPERATING PRESSURE... 1300.0 
TPMAX, Feececcccccccccsvccscvccsccsccce 150.0 
TPNOM, Fe wcs ccc vce cee s oldie os oleic s ollie ees 70.0 
TPMIN, F.0 2 occ 0 0 0 06 0 © © o ole 5 oo sluhstaietetstniale ~631000 
TEMPERATURE SENSITIVITY, it/P. ... 6s es 0.200 
NOMINAL CHAMBER PRESSURE, PSI........ 1030.2 

REQUIRED BURNING RATE, IN/SEC.....cec. 0.615 
ALTITUDE , PD. 2 wos w tie occ ce cece ace eee 10000. 
AMBIENT PRESSURE, PSL. cc ccc cen ccccsccce 10.13 
AMBIENT TEMPERATURE, Rewer ecccrcvcccce ove 483.1 

3 
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TABLE 2 


MICROPEP OUTPUT FOR P, = 250 AND 1030 PSIA: 


Beene. oo Ae, 12% HEEB, 106 Al 


DESIGN 1 Case lof 2 
DH DENS COMPOSITION 
‘AMMONIUM PERCHLORATE (AP) ' ©602 0.07040 1cL 4H 1N 40 
HTPB (SINCLAIR) 130603320 203H.73c 10 
ALUMINUM (PURE CRYSTALINE) 0 0.09760 VNC 
INGREDIENT WEIGHTS (IN ORDER) AND TOTAL WEIGHT (LAST ITEM IN LIST) 
78.0000 12.0000 10.0000 100.0000 ; 
INGREDIENT VOLUME RATIOS ! 
AMMONIUM PERCHLORATE (AP) 70.487% HTPB (SINCLAIR) 22.995% 
ALUMINUM (PURE CRYSTALINE) 6.518% 


THE PROPELLENT DENSITY IS 0.06362 LB/CU-IN OR 1.7610 GM/CC 
NUMBER OF GRAM ATOMS OF EACH ELEMENT PRESENT IN INGREDIENTS 


3.895570 HH 0.878965 C 0.663847 N 2.667428 O 


0.370645 AL 0.663847 CL 
CP = 43.62668610 


ee ee i ee ee 


CHAMBER RESULTS FOLLOW: 


~ T(K) T(F) P(ATM) P(PSI) ENTHALPY ENTROPY CP/CV GAS RT/V MOL WT 
3248. 5388. L701 250.00 -46.80 240.56 1.1990 3.644 4.667 27.44 


DAMPED AND UNDAMPED SPEED OF SOUND= 3211.743 AND 3563.987 FT/SEC 


HEAT CONTENT (298 REF) 1253.397 CAL/GR 2257.631 BTU/LBM 
0.98321 H20 0.71991 CO 0.58380 H2 0.57339 HCl 
0.32835 N2 0.17947 Al203* 0.15900 co2 0.11256 H 
0.07546 HO 0.07476 Cl 0.00880 O 0.00709 NO 
5.93E-03 02 4.51E-03 A10Cl 3.65E-03 AlCl2 3.20E-03 ‘AlCl 
1.32E-04 A1C13 1.28E+04 C12 1.19E-04 AlO .~—-«6.37E-05 Al 
2.31E-05 CHO 1.73E-05 NH 1.61E-05 N 1.21E-05 Ho2 
7.13E-06 NH3 5.55E-06 NHO’ 5.24E-06 Al20 4.85E-06 AlH 
2.03E-06 CNH 1.91E-06 AlHO2 1.09E-06 CH20 
CP = 41.92294693 be 
EXHAUST RESULTS FOLLOW: : , 
T(K) T(F) P(ATM)  P(PSI) ENTHALPY ENTROPY CP/cV GAS RT/V. MOL WT 
2313. 3704. 0.69 10.13 -108.91 240.56 1.1995 3.510 0.196 28.49 


DAMPED AND UNDAMPED SPEED OF SOUND= 2650.051 AND 2952.060 FT/SEC 


HEAT CONTENT (298 REF) 851.166 CAL/GR 1533.130 BIU/LBM 
102958 120 0.67108 CO 0.64916 HCl 0.59218 H2 
0.33181 N2 0.20785 CO2 0.18284 Al203& 0.01785 H 
0.01421 Cl 0.00497 HO 0.00233 A1l203% 0.00021 NO 

t. 338304 ALCl2 1.04E-04 O 1.01E-04 Al10C1 dvoru=es O2 


4.03E-05 Alcl 1.25E=-05 AlC13 9.06E-06 Cl2 
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TABLE 2 CONTINUED) 


Re 
aie 


MICROPEP OUTPUT FOR P, = 250 AND 1030 PS 


Po. = 10.0 BSTA. (83 APY Me eee ee - 


SECOND LINE 


PERFORMANCE: FROZEN ON FIRST LINE, SHIFTING ON 
ISP GAMMA T* _—— P# Ce CF OPT EX D-ISP A*M. EXT 
226.7 1.1835 2975. 9.65 5060.1 1.442 4.31 399.2 0.62924 1976. 
232.5 1.1497 3045. 9.77 5082.8 1.472 4.72 409.4 0.63205 2313. 
BOOST VELOCITIES FOR PROPELLANT DENSITY OF 0.06362 (S.G. OF 1.761) 
5./23451. 10./18585. 15./15857. 25./12612. 30./11520. 55./ 8215. 
60./ 7788. 69./ 7128. 71./ 6998. 88./ 6064. 100./ 5548. 150./ 4113. 
175./3647. ~00./3277. 300./2336. 1000./ 780. 3000./ 269. 5000./ 163. 
EXP. EXIT EXIT EXIT OPTIMUM OPTIMUM VACUUM VACUUM SEA LV SEA LV 
RATIO PRESS PRESS TEMP IMPULSE IMPULSE IMPULS IMPULS IMPULS IMPULS 
ATM SI _k SEC SI SEC SI SEC SI 
1. 9.769 989.6 3044. 104.9 1029. 195.7 1919. 186.4 1846. 
2. 3.256 329.9 2678. 175.0 1716. 235.5 2309. 216.9 2149. 
3. 1.215 123.0 2386. 214.4 2103. 248.3 “2a. 22054. "aimee 
4. 0.833 84.4 2314. 226.6 2223. 257.6 2526. 220.5 2184. 
5. 0.644 65.3 2313. 234.5 2300. 264.5 2593. 218.0 2160. 
6. 0.523 53.0 2311. 240.7 2361. 269.9 2647. <2ieaoT (cane 
7. 0.439 #4.5 2310. 245.8 2410. 274.4 “eso. “Zoges 2074. 
8. 0.378 38.3 2310. 250.1 2453. 27802 Bayes. Zosnomeccnan 
9. 0.331 33.5 2309. 253.8 2489. 281.5 2760. 197.9 1960. 
10. 0.294 29.8 2308. 257.1 2521. 284.4 2789. 191.5 1897. 
11. 0.264 26.8 2307. 260.0 2550. 287.0 2814. 184.8 1831. 
12. 0.240 24.3 2307. 262.6 2575. 289.3 S@a7. 177-9 ime. 
13. 0.219 22.2 2306. 265.0 2599. 291.5 2858. 170.7 1691. 
14. 06.202 20.5 2306. 267.2 2620. 293.6 2878. 163.4) 1619. 
15. 0.187 19.0 2305. 269.2 2640. 295.3 2896. 155.9 1545. 
16. 0.174 +%17.7 2305. 271.1 2658. 297.0 #2812. 148.3 iaaor 
17. 0.163 16.5 2305. 272.8 2675. 298.5 2928. 140.6 1393. 
18. 0.153 15.5 2304. 274.4 2691. 300.0 2942. 132.8 1316. 
19. 0.144 14.6 2304. 276.0 2706. 301.4 "92956. 124.9 1237. 
20. 0.136 13.8 2304. 277.4 £«°.2720. 302.7 2969. 116-9 1158. 
21. 0.129 13.1 2303. 278.8 2734. 30490 2apa. dos.9 lea 
22. 0.123 12.4 2303. 280.1 2746. 305.1 (2062)05100-8 “Saane 
23. 0.117 11.8 2303. _ 2612.3 2758. 306.3 §300anus 92.6 Some 
24. O.112 11.3 2302. "282.5 2770. 307:3 3004. eee Soe 
25. 0.107 10.8 2302. 283.6 2781. 308.3 3024. s7ee1) “ear 
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TABLE 2 (CONTINUED) 


MICROPEP OUTPUT FOR P. = 250 AND 1030 PSIA: 


p. = 10.0 PSIA. 78% AP, 12% HIPB, 10% Al 
S00, ST 


DESIGN 1 Case 2 Of 2 


DH DENS COMPOSITION 


‘ AMMONIUM PERCHLORATE (AP) -602 0.07040 1CL 4H 1N 40 
HTPB (SINCLAIR) | 13 0.03320 103H 73C_ 1o 
ALUMINUM (PURE CRYSTALINE) 0 0.09760 1AL 
INGREDIENT WEIGHTS (IN ORDER) AND TOTAL WEIGHT (LAST ITEM IN LIST) 
78.0000 12.0000 10.0000 ", 100.0000 
‘INGREDIENT VOLUME RATIOS — | 
AMMONIUM PERCHLORATE (AP) .70.487% HTPB (SINCLAIR) s* 22.9958 . 
ALUMINUM (PURE CRYSTALINE) ~ 6.518% , 


THE PROPELLENT DENSITY IS 0.06362 LB/CU-IN OR 1.7610 GM/CC 


NUMBER OF GRAM ATOMS OF EACH ELEMENT PRESENT IN INGREDIENTS — 


3.895570 H 0.878965 C 0.663847 N 2.667428 O 
0.370645 AL 0.663847 CL 
CP = 43.75085831 


CHAMBER RESULTS FOLLOW: 


T(K) -T(F) P(ATM)  P(PSI) ENTHALPY ENTROPY CP/CV GAS’ RT/V MOL WT 
3375. 5616. 70.08: 1030.24 -46.80 230.37 1.1952 3.597 19.485 27.80 


DAMPED AND UNDAMPED SPEED OF SOUND= 3247.693 AND 3603.495 FT/SEC 


HEAT CONTENT (298 REF) 1309.442 CAL/GR 2358.580 BTU/LBM 
| 1.01084 H20 0.71716 CO- 0.59375 HCl 0.57367 H2 

0.32895 N2 0.17941 Al203* 0.16171 CO2 0.07526 H 
0.05727 HO. 0.05233 Cl 0.00586 NO - 9.00506 AlC12 

4.49E-03 0 3.70E-03 AlOCl 3.15E-03 02 2.63E-03 AlCl 

3.14E-04 AlC13 1.83E-04 Cl2...—s-« 7..24E-05 Al0 5.77E-05 CHO 

3.63E-05 Al 2.75E-05 NH 2.71E-05 NH3 1.55E-05 N 
1.33E-05 HO2 9.39E-06 NHO 8.06E-06 CNH 5.08E-06 AlH 

4.45E-06 CH20 3.50E-06 Al2 2.16E-06 AlHO2 


CP = 40.18283081 
EXHAUST RESULTS FOLLOW: 


q 
T(K) T(F) P(ATM)  P(PSI) ENTHALPY ENTROPY CP/CV GAS RT/V MOL WT 
1852. 2874. 0.69 lovigee=i3Gr14 98230552 1.20607 39861 0.197 28.64 


DAMPED AND UNDAMPED SPEED OF SOUND= 2373.392 AND 2644.581 FT/SEC 


HEAT CONTENT (298 REF) 659.260 CAL/GR 1187.465 BTU/LBM 
0.98547 H20 0.66303 HCl 0.63198 CO 0.63022 H2 
0.33192 N2 0.24696 CO2 0.18531 Al203& 0.00097 H 

7.87E-04 Cl 1.10E-04 HO 2.78E-06 AlCl2 2.25E-06 NO 
2.21E-06 A1C13 1.52E-06 NH3 

PERFORMANCE: FROZEN ON FIRST LINE, SHIFTING ON SECOND LINE 


onl 





TABLE 2. (C@MTENUED) 


MICROPEP OUTPUT FOR P, = 250 AND 1030 PSIA: 
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P. = 10.0 PSEA. 78% AP, Is SoRES ye eel 
ISP GAMMA T* px cx CF OPT EX D-ISP A*M. EXT 
- 262.5 1.1907 3081. 39.69 5084.9 1.661 12.28 462.2 0.15344 1610. 
269.3 1.1617 3139. 40.09 5126.9 1.690 13.25 474.2 0.15471 1852. 
BOOST VELOCITIES FOR PROPELLANT DENSITY OF 0.06362 (S.G. OF 1.71) 
5./27165. 10./21528. 15./18368. 25./14609. 30./13344. 55./ 9516. 
60./ 9021. 69./ 8257. 71./ 8106. 88./ 7024. 100./ 6426. 150./ 4764. 
175./4224. 200./3796. 300./2705. 1000./ 904. 3000./ 312. 5000./ 188. 
EXP. EXIT EXIT EXIT OPTIMUM OPTIMUM VACUUM VACUUM SEA LV SEA LV 
RATIO PRESS PRESS TEMP IMPULSE IMPULSE IMPULS IMPULS IMPULS IMPULS 
ATM SI K SEC SI SEC SI SEC SI 
1. 40.059 4057.9 3139. 106.6 1045. 197.7 1939. 195.4 1936. 
2. 13.353 1352.6 2722. 176.8 1734. 237.6 2550. "233.0 gamer 
3. 4.904 496.8 2391. 216.7 2225. 250.1 2453. 243.3 24m. 
4. 3.289 333.2 2270. 2@8e5.° 2250. 259.4 2549. “sora "eae... 
5. 2.435 246.7 2183. 238.3 wW@87. 266.0 2609. 254.7 2523. 
6. 1.913 193.8 2115. 245.0 2402. 271.1 2658. 257.4 2550. 
7. 1.563 158.3 2060. 250.3 2454. 275.1 2698. 259.2 2568. 
8. 1.314 133.1 2014. 254.6 2497. 27815 2731. 260.3 $2579. 
9. 1.129 114.3 1975. 258.3 2532. 281.4 2755. 260.9 2588. 
10. 0.986 99.9 1940. 261.4 2563. 283.8 2783. 261.1 2586. 
ll. 0.873 88.4 1910. 264.2 2590. 286.0 2805. 261.0 2585. 
12. 0.781 79.1 1882. 266.6 2615. 27.9 Zeee. “geo.7 soa. 
13. 0.706 71.5 1858. 268.8 2636. 289.7 2841. 260.1 2577. 
14. 0.643 65.1 1835. 270.8 2655. 291.3 2856. 259.4 2570. 
15. 0.589 59.7 1814. 272.6 2673. 292.7 2870. 258.6 2562. 
16. 0.543 55.0 1795. 274.3 2690. 294.0 2883. 257.6 2552. 
17. 0.503 51.0 1778. =275.8 2705. 295.2 2895. 256.6 2542. 
18. 0.468 47.4 1761. 277.2 £2718. 296.4 2906. 255.5 2531. 
19. 0.438 44.3 1746. 278.5 2731. 287e40u2omee esac Dana. 
20. O.412 41.6 1731. 279.8 2743. 288eqmecceme mess. 0 2506, 
21. 0.386 39.1 1717. 280.9 2755. 299.4 2936. 251.6 2493. 
22. 0.365 36.9 1704. 282.0 2766. 300.3 2944. 250.2 2479. 
23. 0.345 34.9 1692. 283.0 2776. 301.1‘°2952. 248.8 2465. 
24. 0.327 °33.1 1681. 284.0. 2785. 301.5 2960. 24753 2a50. 
25. 0.311 31.5 1670. 284.9 2794. 302.6 2968. 245.8 2435. 
+ 


The output from "MICROPEP" indicated the following 


* 


C = 5126.9 ft/sec 
th 


I = 269.3 lbf sec/lbm 
SPth 


em. 13225 
Moles of condensed Al 503/100 gm = 0.1794 
Pp = -0636 1lbm/in? 

The AGARD~-AR-230 empirical equations were then used to 
estimate the value for nc - Nex was chosen to be 0.93. For 
£f the calculated residence time of 7.2 msec, Figure 3, indi-~ 
cates that a better value may have been approximately 0.87. 
As discussed above, the losses are iterated with the throat 
area. The resulting losses and throat area are presented in 
Table 3. 


Once the losses have been estimated Isp can be deter- 


mined from 


Isp = cae ea 


The following parameters were then calculated: 


m EV ANe 

Vo = ee eee 

Ve = VD/Va 

Doe = (YI-V_)D 

Ap = (Z)Dpe = FD (1-Vy) 
ee eh, an 


These parameters are shown in Table 4. 


mye 


PERCENT 
PERCENT 
PERCENT 
PERCENT 
PERCENT 
PERCENT 


Etacf= 

PERCENT 
PERCENT 
PERCENT 
PERCENT 
PERCENT 
PERCENT 


ESTIMATED PERFORMANCE LOSSES 


DIVERGENCE 


0.924 


TABLE. 3 


LOSS is 60.6 6 6 6s 6 ote amenenener 
TWO PHASE FLOW LOSS... ccecccccces 
BOUNDARY LAYER LOSS... cece ccceee 
KINETEC BOSS cc 2 oe 0 0 0 « 6 os eleienelo sie atene 
SUBMERGENCE LOSS... ccc cccvcccccece 
EROSION LOSS. cc cc ccccccccccvcccccs 
THRUST COMPFICIENT EPPFICBENCY. . < 20ers 
1.562 Ath= 
DIVERGENCE LOSS... cece rece ccc ccvce 
TWO PHASE FLOW LOSS... eee. ee cece 
BOUNDARY LAYER LOSS.......- eee 
KENETEC LOSS 2. 2 2c 0 co 4 sete ce Mieletelste slate 
SUBMERGENCE LOSS... ccc cw cercrccccvcs 
EROSION LOSS... cc wc cre ceeccccvcces 
THRUST COEFFICIENT BERPICERNGCY .. .2 . setts oie 
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1.704 
3.202 
1.828 
0.849 
0.000 
0.007 
0.924 
8.856 
1.704 
3.199 
1.783 
0.849 
0.000 
0.007 
0.925 


: 
| 
| 
| 
: 





TABLE 4 


CALCULATED MOTOR DATA 


THRUST COEFFICIENT... ccc ccccccccvvvccves 1.563 
BUI eo meee wee see eee ease ee sce sere crees O29 25 
SPECIFIC IMPULSE, LBF SEC/LBM......... 231.587 
MASS FLOW RATE, LBM/SEC. oo cc eevee wee 61.532 
PROPELLANT VOLUME, CU IN... eee eee ween 6286.8 
NOZZLE THROAT AREA, SQ IN... eee eeeevees 86851 - 
NOZZLE EXIT AREA, SQ INw~- ccc ecvcces Liye 0 
EQUIVALENT PORT DIAMETER, IN... eee. 7.476 
PORT AREA, SQ IN cc cece cece cccccccvcees 43.893 
PERIMETER FACTOR... cccccccccccccccccces 2.489 
SMMCET) WEE, ENw ccc tcc ec ce ec eee sce e ee. 2.500 
REQUIRED WEB FRACTION. .cccccccsvvevaces 0.312 
VOLUMETRIC LOADING... .ccccccccccveccces 0.782 
LENGTH TO DIAMETER RATIO... ccccccccccce 2.500 
PeeeNGE Permit, MSEC... cc ccccccecccvccs 7.2256 
THROAT TO PORT AREA RATIO (J FACTOR)... 0.202 


REQUIRED BURNING RATE, TN) SEC... 756 5 o 0.385 


2 


At this point the following expressions remain to be 


Satisfied: 
Perimeter Factor = PFAC = Pwetted/ "pe 


PFAC =Ap/™Dpeb 


Ap = Pwetted 


web/*, 


cnom 


E 
In these three expressions the unspecified variables are 
PFAC, Ap, ©v, and web. 
For a selected web or r the other parameters can be 
determined. For example, in the current problem Table E 


summarizes the required values ies selected web 


thicknesses. 


TABLE 5 


REQUIRED PARAMETERS FOR SELECTED WEBS 


web(in) r_(in/sec) Ay (in?) PFAC 
3.5 .538 1796 1 soit 
3.0 461 2095 262 
2.5 2385 2514 2.67 


Once the required burning rate (rc) is determined, appro- 


priate values of a and n are required for the 
° n 
expression r = Po. - 


The grain design problem is to find a geometry which 
will satisfy Table 5, together with the required conditions 


on L, D, Port and Fre Candesemee-)). 
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In the present program, only radially spoked star grains 
can be analyzed. It may not be possible to satisfy the 
motor requirements with this type of grain. The process 1s 
essentially one of trial-and-error. 

Three attempts at meeting the design requirements are 
presented. One 6-point and two 8-point designs were tried. 

A typical output from the "Design Mode" of "GRAINS" 1s 
Shown in Table 6. 

Once the "GRAINS Design Mode" has been run, "GRAINS 
Burnback" is run to obtain A, versus web. 


The output from "Burnback" for Case II is shown in Table 


iz 

"ROCKET" can then be run to obtain the pressure-~time and 
thrust-time profiles. Input data for Case II are shown in 
Table 8. 


A small section of the output from "ROCKET" is shown in 
Table 9 for Case II. 

A summary of the results for the three cases 1S pre- 
sented in Table 10 and pressure-time traces are shown in 
Figure 7. 

None of the cases attempted met the desired design 
values with sufficient accuracy. Further iterations may re- 
sult in a better solution. However, the star design has 
considerable sliver with the resulting long tailoff. This 
section was presented in order to demonstrate the process 


required to obtain a final grain design. 


i, 


OUTPUT FROM 


Mayiague Number of Sookes = @ 


Cesired Port Area = 43, 893293 
RI} = 3.50006 Se ia 
SRi1 = 0.25000 DEL: = 
SRI2 = 0.25000 DEL4 = 


rINWEB = 4.00000 WESSTP = 
CiRTOL = 3.10000 °TioLl = 
Tolerance Ligits 


TABLE 


6 


brain Quter Nadius = 
Desired Perimeter Factor 
2.00000 
000 


0.25000 
0. 10000 
0.17000 
Q. 29000 
0, 10000 


FIN2 
FINI 
FIN4 
KSPOKE 
Per 


Foye Area: Miniaus = 39.90364 - Maxiaus = “5 
Perimeter: Minimum = 92.6139: Maxiaua = 64 
Circle oo: Minimua = 324.00000 Maxl@ua = 336 
med fasckness = 2.5990 
n APoR? tenis rFat RRL Re R3 Re 
CRRUR -> Gi Greater tnan PITCL, theretore edual 
31 682 Sond Oe ee ee 
8 42,379 Sore de Sle ocd cee eerie 
@ 42. Coweta 2.57: calls Sica ee ees 
cance ote og, Of eater san 7 1b, therefore edua: 
So ofleatl Sageee co5nd cilcy  o-0G) MeO. Gy duleS 
seeQR -> €12 create7 tnan -lidi, therefore equa: 
6 22.52: $2.602) 5.143 2,500 eo. /aC eeeae decccmee 
cRRUx 7? ua Greater than rid, therefore egua: 
cRROR <> ec greater than PESTaL, thereiore ecual 
B 2/65! 63.043 2igse. 25209 2ad0 ou oacde mencOrmme 
ceeOr -v Fai4:£ O.02 tess tnan 0.0 
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TABLE 7 


GRAINS "BURNBACK" OUTPUT FOR CASE II 


Soren outs 
maxigco Nugoer Gi Spcaes = 49 arain Juter xacius = ~~ §.60000 
vesired Fort Area =  4$2.22200 Desived Perimezer Factor = : 
Rii = 5. 00009 del2 = J. vvs00 ring = 0.06600 
SRiL = UV. 2d000 Dei: = 0.00000 rikl = nae 
sai? = ce 20008 Doi + = 0.00000 riks = ve 2ul0G 
rikwep = 0.00009 woBSTF = 0.00090 KSPCKE = 0.00000 
CizichL = 9.00000 PriCiL = ¥.v0009 PERTOL = 0.00090 


brain Lengin = 40.00000 Pig = | 2oecueoc Fal = <S.c00 
Sees) os Sie2 20= Si erse= ee aoe ka = ae. 
Number of Spokes = 8. 
wea Thickness = 2. 00000 
Wed Perimeter Perimeter Burn Port 
kurned Factor Area Area 
(IN) (IN) (SQ°TN)  (Séiei:.. 
Burnback of Seki (Sooke Corner Radius) 
9.000 38.737 gecue | ONG! Shee) 
7.900 52.001 2.48: 2962.935 45.068 
0,100" . 33.eba 2e4i5 2374.604 46.027 
C.ia0° ” goers wead? “ceaismmo —Deecds 
R04 (mula. 299Iaias OS) doeoes zs : 


a cE 


J0rd some tau 20905 2585,997 84.747 
ves ee, ie 0% GcewUcugos Nuno dic 
Lace eee 10bBS coUbee/ 2 Oimuas 
seoce Wo eraes acUsc  cd0Crey |) i eembos 
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TABLE 7 (CONTINUED) 


GRAINS "“BURNBACK" OUTPUT FOR CASE II 


ee et te rs me 


1465 $9,117 1.447 2364.680 132.617 
1540 £70426 3,386 2296.950 136,593 
-.613 £5.73 1.328 2225.220 140.247 
Mieeh 53.957 . “eer, 2251.489 142.777 


~ 


Geo 865%. 037 ye 26) 459 162 
348 52.62 10885 208:.855 152,929 
PeCit “52.563 sures 62053.309 «= anges 
a! c/566 1,419 2089.449 165.195 
$5337 36—s- $330 1.99% 30677.196 177.695 
ceauu ce O45 eecee 257.3) ete les 


Percent of fuel remaining = 4.40 : 


Ben iol 


Burnaacy of renaining fuei 


Pees2 820 ROS Gi Fee 26. 073 
Bese) ere OG GFN 155005: 
3.066 6.419 6.208 412.753 199,697 
guc00 Sere STS eh <07 «apo een 
3.343 9.000 96.000 «=, Ge 28 ass 


TABLE 8 


INPUT DATA FOR "ROCKET", CASE II 


gr rrr 
ter Pertorgance ©7dardak [PMP > version 1.6% 
INP SInK oaeene 
MOThS KAS ¢ GRAINS 
DESCRIPTION ch GRAIN 
BURN SATE cr. o946 


BURN RATE EXP. | 9.4000 
Pi SUB K §, 2099 


URN Rer. PRESS. i030.2 
C STAR 4760.4 
: PSAP cht. J, C061 
DENSI VY ), G69 
IGHITION TIME #,0000 
BOLTA Ich. rahe Ue vend 


Norgai inout is finisned 
General Conticurazion PARAMETERS are ..«. 


Tarcat AREA Secu 
Exit AREA fi7e29Cs 
Expansion “4Siic pet sais 
Total Motor VOLUME ee 
AMSIEN, Tepoerature 7Qet 
AME: EN, "ressure 


o- « 
= 


o-= CJ 
cm on cS 
= 
eee 


Closure ELowGu; 


Pere 0 
Landda 2 0090 
Ms C3246 
3afka pee ue 
throat DESioh cressur cee 
Theda: E805. bELAY TI 2000. 008 
CSSSsupe VS. THROAT EROSION RArc 


° 600, 0.000000 os 


we té 


WES BURK AREA 
0.000 2349.472 
0.250 2411,952 
1.108 7507.744 
1.540 | 2236.212 
2.011 2060. 360 
2.500 2097. 107 
3 2,832 812.473 
| 3.342 0.000 


| TOTAL PROPELLANT WEIGHTS 4 IGNITER WEIGHT 0.00 


ww 
cr: 


a 


> 
is 
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CASE II 


: BROCKE. 5 


TABLE. 9 
SAMPLE- OUTPUT -FROM 


ite 


1 + \ a ee 


<> = <> @& i > Coc C5 5h J <)> @W MH so Oo So De oo Oo Do oOo Oo TOTO OlO CO oO oO oO SO OST UD ltt 
=o Co © @2@ @ so Cs @ MH 13 SS GD SO <o SS SS 6 ts? SS USF! Oe hUwShlhOGmhUmLCllhlUG! lh SS a eae ee OC ll a a | a OS Ss es 
= &2} CT © 2 2 t& 2 49 3 Ws > 2 2 & > SS C3 SS aa a ee eS ee eee 
2&2] © gaa] 2 2 OO Oo 82 DS © Ce fo <> 2 45 GD SS 8a cp CD RG ae ee re eee eee es SS 
oe e s s a a e e es a e Ld] = e @ s e es eS = es a e e = = = s a 2 e e a e = co es s es e ™ = 
= © © © © <> €o 65 20s 25 <3 ooo es 6 SS eS SS SR Se OR SE sp st ee ee eee el ee SS ED Cen eee a es Ce a ees 
ot re es lw Olle te Oe em ee we ee ee He ee eee lee ee eH eK SG He we re Uw Hl lll HE lll ll lel r4 ta ae ee ae ae ee te, Si, cee cl 
v7 in Win bo WI wf wD ie iT dso) ie ie) a OU I Ww Wl Won oo bo 
onmrqaouounenocnmnn®s - wm wrocd rn «(oO mM ro wo wo Mmm ww wo mw mw wen co 8D co tn os www Ow Ow wD Ww co loo lo Cl 
= es @ e e G e e e = s = e e es e e e BS e ea a 3 = es es s es a e es e e e es 2 s 2 e = e e 
mo wo wmcrmocmimnmn mm co my (nN 197 CO tN tN 6s LO) oN co ON fro om ws ooo co to lUrtD lhuco lUvcOoOlUrmwmlch!lhLUuem!lLCUue oll Clo ltl ClO 
- 
9 U9 (DO 69 OO 100 8 OS UD Pe ore ct Oo Om. OS Oo) SS to eee Owe ee De ee ao ofr «tr Od OC 64 ON 69 O84 OD ON Om 
90 UD ™ ON 1D CO 4 BP hs ome FF OH 8-8 Om) 0nd nee er IO Om Oo SG Ore ™ SO 9 Dm ogg a DP rd 
e @e e = e ee @ e @ =z e e Ss a e @ 2 e es @ e a @ @ oe e e e e e e e e e e @ BS e es es a e 
wr OD Ora ON et ID OD ee eT a 4 TY CHD bee ed UD tr sk Od od ee ee OOD SS Mw wd ™ 1d) UI IM sr sr 4 
C7 ¢9'¢9 €9 69 €9 C4 OO’ ae Ce © 4 Og 404 ¢-4 €~-§ C4 ih Ce we we a a wee tle leh Tle lO «-4@ 4 ~~ 
3 Oo Oo 2&2 OO CO F232 TW OC DOD 7M OO DW DW OW Oo OU OO FP DO OS OP Ee OD DO DO FwWD CO Oo Oo COCO Oo OS OS 
ao Oo Oo Oo CO OO CO CO CO CO OO OM DS OD cD ed > CO Ooo Ds Mw OF OW YM DW OD Dm DOOM CO CO COCO ff oO Df Oo aS 
e e = se = 2 So e = e e e B es - @ = e e e e es es e = e - s s = e e e e e e s = s eo = s 
oOo Go vo ovo CO 2 Oo Oe TW DW FH Dw OS © DP L “DS cs £9 > SP 6? DS «2 fo SS SSS SS Se SS ES Se a ee 
\ 
wt M™ CA OO WM UD th SO ese AD Ne CD 8 eee OO we OO OR OO wt Ot say ‘Mime TE OH Ow a + we OPH 
OomM Wd rere © CD OH ON CF AN OD ee DD 03d CO PU aD Pr FAO SD ON er Os OS SOU eee eH OT Oia 
e = e = 2 @ 2 a ‘* @ e S e = = = s e = e e 2 = s e s s = - e - @ e e = = e . @ @ @ = e 
CD 4a wee CD 0D OOD Ow Pm 0D UD eR Or OD de ee DEM CD ts Ot TD Da a ad wD ee a es me a do cr 9 
co 09 OD 09 0D Ba nd Oa Oa Cra Oa a rd Crd Ca Crd Oa wrt we 4 et we me be 4 eee oe wt ut ay 
f, A 
wom © wes re Fe UD 4 OS OS Rr te OD ee ON OM O™| os OI |S! un wo ad ome OI SOON mes OT tc ako le 
mw O10 C1 OO we |S 1d Ode ON OO Lgl tr) oes Or ON ad swe et es OD 99 os ok las 
es s 2 2 e s e e e e e = = e e e e 2 es = e @ & e e e s e = s 2 = e e - e = e s s e s 
72 <> we CY 049 09 fh 8k UD FD 8 tO te Pm HD COD Te le wee 8d £2 OD er od I i on oes rE! POU oD aN ee ee OH OO oer ae 
womnmmmaiui@mm awa cow tN 9 oD CO CO cre eT ON Ce Crataty CY CM are rm oy at ere tin Om Am OD CO SG OC CO fC mS D> 
eo 69 09 69 9 OF OD Os OF O79 O93 OD Ulm td Oo le OD OF Oe) eo) oe ew oD i OO 
Come Dea @ HE O&O we 4 9 10 TN 6 ao 3 ee ON 4 978 9 10 Oo Oss PR re ed 6 SS OOS ~~ OG CO DWE Om Oo 
ow wo rf WO CO F 1900 SF OO UI ers tr Om SS a sO om ea ae Oo est i OO 
e a 2 @ - ee ae © e e e e = se e e e ® . 2 @ e e s e e e = % e e e e e s 2 e = s rt = se e 
CF 3F ek SR KK HP CD ON ee DD OO 1D 8k AD Ps Oa a we 0D ot Om ew ed OO et OD wk OD ee Re oe 8 OD OK 
Trim worm OO ON OD wwe 8 OD OD Ph Oe CON Owe OD Or ese Ce Sl wt OD OD sr de owt Er OO Oecd 
ud 69 WI 1 wD OD 8 kat 11 ta td ae eos Re Oe PS Le ese Pe ss se oe 0 oo cod oo tt see 
P~ Ps Pe. Fe Pw Fee Pew Po~ I™ Pen Ps fs UF T~- Pes few Pew Pre Pn Pee Pow Pos f~ F~ P= f<. Fs. t= a es Pa er = ~~ eS i ES Or Pe 
‘ 
woo IM es OID? INW ON ME SO wb) oS Oe co e1 OO CHA WY Han ceo wm wo OW OF Ys HINO MIN Ow CO UD WO 
“TP se om OIF OF 1A Oo EN od ee enn mo Nt HOT Oe OO in OOO Ol aN Is. Ulett Uh row oll wr Or lO OUT ee Uc OS rt 
e e eo = © @ € e e e & e © e e 8 e e 2 a e e e e -e e e . e e . e e * e @ e e e e e s 
Ud wr or CF Ce we COO EN ££ OD 18d Pre Pee OU OD gt > 6-8 FD ee OD ED OF we Pe OD SF se DD Ue CO OO er WO OO tr oD ON Ss CD 
sF OD Oa wee OO OM CO 8H Ud Nr OD Oe we > EY At we rd od ae OO a i 80 8 in aoe 
Wao woo 3O DD 3 OD KR KP Nr OO) ON OD ed OD Hd ad Od Oe Hd ON Of f- 46 t § $ 5 $ t 
$ f 6 f o ‘ § t 6 i 6 I q 6 #” $ ' t $ t $ 6 i ' f 1 $ 6 6 i i 6 
uscd ud Od wt oes OT OE oar 0 1 wo et a ad oe te te Oe OO mM om i 
ao re OOF -~ FP 0D UD F™ OO 10D PR DD Oe me TO 8 i od ee ON ON ON NY INO OD oOo CUCoOlUmmlhLre TD let OO 
e e 2S = e @ e 2 e 2 2 e os e = e e e e es 2. e e e e @ e a e e € e . 2 e e 2 e e a = e 
CY ON a Oe Os OD OC 2D cy 8y Om 09d CD 2d UD PO Oe er OE CO Uo ae tN OO ei wool) Ocoee se re OO OW 
tr IW er To Cd SS OD r~ 1D ws reediwe SS Ce rs tO rd me OO TY Ps A OSE OOD UL oe a ao ee oI UNGS Oa a 8 Cl 
u79 I wr er oe Ar ek OD 9 OD OD OD OD OE ON 8 Oe OM OO wre wet wt wt eet et we we 
w9 «wo AwO fh Fe om WTO MAD DW le ote ee MN OA OD er UD oe wo NHN Oo ee OE OI et OD OOPS OOM SH TT OOO OOOO OO 
ctw OM SO CA wr 10 OC MH) bt te Om ee OD Pe et OO Oe ee OO OE OG oD ODS COTE Or OOO Da GS OO 
cw OE ON ed 09 9 9 0D Rd ee ek UD 9D ID toe eee er OUD Ol os swt we ONO Oe Pe ie 09 OD so 
2. e = e e @ e e 2 = e es s s es e 2 e es e s a @ = td 2 2 e Bs = a e = = e @ e s = @ es e 
~~ Fee Ps Re Fe Pe RS Pe Pee Roe RR hs se PS Pee Pf. Pe Pr RR PR rm Re Pf. Ph TF. 1 wom wowcaoreo wo om aw er ct OY on 


9.00 
9.00 
9.00 


2.64 
de 
et 
aarti 


@ 6 


4 
e 


403.50 
$14.6 

Pod 

63 


CAA £ 
ova7 ob 


See l 
eG 


“74 


$4.65 
48 
24.48 


10.062 
10.383 
ey ag 


5 


AN, 7 


L 


2 


Bo oget? im 
2, 


~ 
= 


A 
ae 


a 
a) 


1.300 


* 


Os 


CASE ITI 
daeecd 


A 
2 


INT F 
86425.01 


“ROCKET y, 


wo. 
a 


am = 
ae 
dewe 
¢ 


BV eee oO 


17 


iNT PC 
5906.01 


vu aes 


TABLE 9 (CONTINUED) 


4368.49 


_p~r 


NARS 


SAMPLE OUTPUT FROM 


“a 
we 


4 


BS os secs 
Cowen 


rua 


é 


AVerAnz P 
MAX 

Fe 
379.55 


iahe 
7.573 


we 


poe cme ay 2 fr. PY FX 


im ae 


~~ i" mica? 0 Ce oD Civ ed UES D KD CD OV AD : i 
SO OS Ob Oe ea 8 SIS oS Ww 8 re ID FO 10 i ee a 
a B e s s a s oe eo = e 

x, ™ 5 2 4 = « See ee 7 . 7 S 3 re 5, oe es = VS eee | Ace esi a ry at ACN = ss ri = 
oe " es - io . jeep 2k oe) €-8 O88  w 8 C-2 F-h s e~d 2 8 we ar ~ ° 
* Ae cy = 4 : 40-8 Ca ad’ ate Oo C4 CW Ea Ce a TE a OE Oe Oa Oa Od oN 
in~ > OF ow a ca ee Ca me 8 CTD oot 8 WY NR Cam oe Oe at td tee KI > €Mm ob rm 649 
Od UD OS le Ree RO Se bed ome Fe Oe Zt we bea e le 2 ek 8m 8 8 a OE SO OO 

oe os eo *. 2 a = ea eo e es e = B s = eo s se e s e = es oe e 
Z> we C2 OD C8 Fo oaea ded eee Oe we Pom 0d ee tt ee ea fe § a oad Fe OTD AF 
By CRY Che ary Use Ke Ay ae Fe nee CY CN PY CDOT Zy ep SS ete NLS pes coo tag) ee | COC | Ge 
oI C9 C3 OD 0 Os OD C9 CD 0 OD Oe’ OD ey eg wf 6 F st °F 32> 4, «6S COKE: 

= 

Pr. Fee aap Sa or OD et 8 8 OO ON Oe Or a to OR Oe eM Oa oA 
P-~ 0c3 OD ef CD wee cm tt CO rs 1 Om bd ow ti OD OD we Pe cf Cry ON OD 
Ce CH a OM Ca Ta’ Ot wt ee me we oe te ee DD CD UC EN Os OD ofa f- 

® Ss es e e a e a e - e @ es es e s s e e Ss e e = es a e 
w¢ ~ «4 ¢ ~~ 4 -~@ @ 020 @ w4 . @ « @ ot 4 r @ w 4G « «4 a | _—4 1 4 v ¢ Sar CS = —s ai a 
4 wD a 8 OO 1D 008 AOD Pe UD Me 0a em 2d OS OD 0D ced 1 be ee SD VP Tr 
sree DT TD dD os oe 9 oe 

es a eo eo s a Ss oS e eo es > es ea es a @e sa o es a BS Ss @e a eo 
--~ F808 C9 8 “3 007 8D a Fs P™ OD FO Rem oD ewe Oe OD ot od cad FR OD ots oh 
Cty Ory CTH te ay he ate EF EY FFD LY oN CY CPT «D> st RES te hI KS ce ll SCS CS 
m9m— a ss I OD FD ed DOD Kee <ge PF «xe OSS OA Or oe <3 st °F «fF 

= 
‘ gf 

m'co wn 42> 60> PS CO Ff SD Co SD Wm C9 UD te th we a SD Owe ee ITD ON 
we oe wr Fr LCF LT CO wt Pe ee Poe PE OD 8a 8s SD CD LU ud nd OE Oe UO Owes 

oe a eo @ a s gg es e @ e e ee = eo eB es oe os es S Q s s B ea 
CO OM me OD 64 ED PS ED a bree OD tte SD FO ery add od ee re OD Ar ORF Wee 
coo er ev) C4 WI €- OO CN CO £O oO Sr 078 Mery OD re I ops 0) fen ee =r Un PR. 1 
ud &. CO C> wee C4 OF hr OD 1 eh ON ON mete CY CN oo Ps. oem oct Cd 09 SP 8S 81o8 
wO 160 10td F™ Rs OT ie Ee Bee Pes ~~ OPS ORs Pe-0O 09 ots Cn) fn) 0 ty fry Cty Cm ot ety 
cc 9 (9 OO ADO OD OF) ore OF oe CH fe oO kd TO Ur oo ate 00 «69d OC UD hum 
OOo wd 83 dD DD aH De AO Cited ee Dg or OO os > co (Dp «xt 
wd oe O48 OD OS) CH ta OD OO A eee a oe Oe oe oe gee PR I Om ots 

es es s a J e a oe e 2 oe es e e s 2 s a a e a es a 2 a e 
UD em Oe we aD Fs Fe Re RO ter oa 0 ade et PR. oo ore. YP eee a co od tr ON 
CD 82 a Os. OO OD ee ad OD ee oD UD ep. cry et oa’ ee AO ORS. ED eed C9 
wy is wd DD wD I Ww 1A 1D 8 1 td Oe Oe OR OS OR CO UO Un CO 
9 asd dt a OA i wb 8 nt OO tthe 1 ee 8 tt 8 Oo lwp kn to 
om er oD wr re. s 209 602 Cea Bs tte sk Re Oren ee ad SD Od SD AO CD CO CUP} om SS 
co 4 OO tr fr 0D Dee CQ TD OD 88 wee w-e OM OD Or Oe ONO wee 

@ eo Be oe oe s es es es e es e e e es es a @ es @ s es * a a e 
CoP oad IV we OD 8k SR eC ee DD OY Oe Ca es ea OD Oe Oo doer 
Pew oF © $F CF CM CD OD CO oe +e: (4 CO Fy. Ced CS UD wee COV ed LO be OD Che ONS 
A aD oh C4 Sl CUD CK NP OD see eS rem. 8s 8k od ee D> oe 00 tr OD Od 
Con CF CD CF OD AO Oa OE Ca Od Ca ee we uate pee emt le Oe 


‘owe OD COP 8 OOK SD NE OO Od te wet eee OD OO ee Om LO Owe tp eee fA OOD 

Co f°) 6ST UMP ONS o7Fe7 wow Oo re~ r™ CO TDN ay @ «4 ) 82 69 DD WO OD OGD ON 
es os * oe s e eB a es oe es es eo ee s @ a eo e e » @ @ @ ea es ® 
oO 69 OF COD CF OO OF OHO OD OY FD ed OD Ey CO WR Or er RR Ot Oe 8k oe OD od 
eo TH Cs > 69 OO WD Pew Few de OD wee ee ae a rm rm a tO on 6 of 
10 C4 OD sad 69) 6 te wed YN TON Cat? D> Ch wmwmeoi nmr CF +e OF 2 tt 
es a e ea s = a s e e es es @ e e os a es e a ea e e es 2 a 
7-2 te DP CTD coe Fe 4 oF 2M OY 1D 160 60 a CoD oft ome rs 1 tO Wt or ot 
“ae CO Ce IY OD or OD kD? OO A Os 8a OP Ot Nw me me to oes Uo lle Oo 

Ce Cxg os C-4 Cf ° ve 


2°@e-4 et wt wd wre 


rT co oF °F UD OO ~~ COC ~-<« ~H C4 ww or et ch WD 
wom CD RP ee OO DO YF 1 OM 6 DO Or Ar LO md oO? wt OO} ON Or 
eae MVD SP we ee ee OD Se to. Ure we C2 UD OR 

BS e 2 e e 2 @ e e es es e es 
r~ Pe. ot. Pe Ps F299) OO LO Om Curd CoD nm toy CA TY are re tS 


0 
16.389 





0.345 2l 216.93 


411.20 


85779.98 
64 


6875.28 


5056 
43.9! 


6.13 


24.48 
15.62 


liooin 


TABLE 10 


SUMMARY OF RESULTS FOR THREE GRAIN DESIGNS 


pa — 
GRAIN 3 
No. of pts. 
web, in 
Required L (in/sec) 
Polk. 
Port Area, in? 


RR1, in 


Churn » sec 


total 


16 IN, 


Crs 
Die 
42.48 
2.750 
Zero 


LOmOs 
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12,464. 
6,956. 


Ores. 
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FIGURE. 


PRESSURE-TIME PROFILES FROM “ROCKET” 
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X.- CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

An interactive computer code for the preliminary design 
of solid rocket propellant rocket motors ("SPRMD") was suc- 
cessfully developed and its use was demonstrated through a 
design example. The code incorporated several existing PC 
codes and AGARD recommended empirical expressions for losses 


in aluminized propellant systems. 


Be. RECOMMENDATIONS 

Further modifications to improve the capabilities of the 
code are essential in order to make "SPRMD" a more effective 
preliminary design tool. 

- There is only one grain configuration (radially spoked 

star grains) incorporated into "SPRMD." 
It 1S necessary to add other options for grain design 
in order to increase the grain design flexibility. One 
possibility would be to incorporate the Grain Design 
and Internal Ballistics (GDB) module of Solid Per- 
formance Prediction Code (SPP) (Ref. 3). 

- "SPRMD" also does not currently incorporate any plot- 
ting capabilities. There is definitely a need for this 
program to have plotting capabilities for the grain 
design results and the thrust and _ pressure versus 
time profiles. 
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"FLYIT" was not incorporated into the program, although 
"SPRMD" was structured to have flight simulation capa- 
bilities. It is also recommended that "FLYIT" be 
incorporated into "SPRMD" as another enhancement. 

"SPRMD" iS a primative type of preliminary solid pro- 
pellant rocket motor design code. Further enhance- 
ments, such as propellant searching, screen editing, 
and increased solution speed for "MICROPEP" and 


"ROCKET" would make this code more user friendly. 
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program SPRMD.- 
CRRAKRKRKRRRKRRKRKRKRKKRKRKRKKKRKKRRKKRKRKRKRRKRKKRKRKRKRKRRKRKRKRKRKRRKRRRKRRKRRRRRRRREK 


* 
SOLID PROPELLANT ROCKET MOTOR DESIGN * 
THIS PROGRAM (SPRMD) SOLVE THE PRELIMINARY DESIGN PROBLEM* 
OF SOLID PROPELLANT ROCKET MOTOR. SUBROUTINES INCLUDES: *# 
(1) INPUTO: BASIC INPUT FOR DESIGN * 
(2) INPUT1: CHEMICAL COMPOSITION INPUT * 
(3) MICROPEP: PERFORME THE PROPELLANT EVALUTION *. 
(4) PFE: SIZE THE THROAT OF NOZZLE PORT AREA OF GRAIN ETC# 
(5) LOSCF: CALCULATE THE THRUST COEFFICIENT * 
(6) INPUT2: STAR GRAIN INPUT * 
(7) GRAINS: PERFORM THE SPOKED STAR GRAINS DESIGN * 
(8) RCKTIN: INPUT THE DATA FOR THRUST AND PRESSURE PRO- *'° 
FILE CALCULATION * 
(9) ROCKET: COMPUTED THE THRUST AND PRESSURE PROFILE * 
kkkkkkkkkkkkkkkkkkkkreerkrkkkkkekrekkkekekekkkkkkkkkehkkhkkhkhkaaaka 


gqagqgqagqaagaagaagaadadgnqaaanaa 


real favg,tb,1,d,meop,pc,rrqd,pamb,tamb,h,ld,webrqd,wf,m 
common /micrp/ a(12,12) ,kr(20),amat(10,12) ,jat(12),in,is, 
1fie(10,6),ie(10,6) ,alp(12) ,w27,n,dah(10) ,rho(10), 
2wate (10) ,wl(6) ,w43,ig,np,vnt(201) ,w47,name,ser, 
3floor, itag(100) ,wing(10) 
S integer 
character*1l0 yes 
CRRRRRRRRRRRRRKRKRKRKRRRRRRRRKRRKRRKRRRRRKRRRRRKRKRRRRRRRRRRKRKR RRR 


c * 
c BASIC INPUT FOR MOTOR DESIGN * 
c : ® 


CREME KKK RK RRR RRR 
1000 call input0(favg,tb,1,d,meop,pc,rrqd,pamb,tamb,h,fmin,tpmax,tpnon, 
*tpmin, pik) 
write(*,1)favg 
Leen , PvErGe THRUST, LBP. cccscvcecsscsccccese ',L legen) 
write(*,2)tb 
2 format(lx,'BURN' TIME, Cie OP scutes. «ccf oncksiyee ‘yp flv2% 2 ) 
write(*,12) fmin 
LZeeeOGnee@ieaMIintMUM THRUST. REOULRED, LBF. « sities oom ‘yedi2.6) 
write(*,3)1 
See OPemeeeiee ee er LENGi ,  LNCHicssie.c « claisiene coccoccess ',Lide2) 
write(*,4)a : 
4 format(1lx, ‘MOTOR DIAMETER, mame mts tare. 2) 
write(*,5)meop 
5 format({lx, 'MAXIMUM EXPECTED OPERATING PRESSURE... m2. 1) 


write(*,13)tpmax ‘ . 
3 format (lx, 'TPMAX, ee ooo sueMeinialMsMeienetcii isle alec © oi sc ele ielele sate. c « 1) 
write(*,131)tpnom x 
Pw OGM ieec Tere, Fence cco e esses esse cecscccscccsce ',4241) 
write(*,132)tpmin 
TS 2 pO Ct, DM, Fes sestcccccevecsscccosecsecscscsse ‘',f12.1) 
write(*,14)pik . 
14 format(ix, 'TEMPERATURE SENSITIVITY, @/F.ceccceeee *,€12.3) 
write(*,6)pc 
6 format(1x, 'NOMINAL CHAMBER PRESSURE, PSI........ ',f12.1) 
write(*,7)rrqd 
7 Ef£0rMAr ix, ReeerneD BURNING RATE, IN/SEC.-cocceee *,B12e3) 
write(*,8)h 


g format(ix, "ALTITUDE, Palas ®eeseeseeesees1eeessese#e#seeseee#esss#seese#s8se#?*? : # fie? @ 0) 
write (*,9)pamb 
3 format ( 1x 4 ‘ANBIENT PRESSURE, PSI @®eeesee@seeeseeee#seses8ess8ee0ee8kte @ ‘ # f*L2 e Pa ) 
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write(*,11)tamb 
ll format(1x, ‘AMBIENT TEMPERATURE, Roc. « s oe seisteivisict pier) 

write(*,10) 

10 format(1lx, 'ANY CORRECTION OF BASIC INPUT? Y-YES, N=NO... ') 
read(*,20)yes 

20 format(al0) 
if(yes.eq.'y'.or.yes.eq.'Y¥') go to 1000 

e 

CRRA EAKREEEKEREKKEKKEEER 

c * 

c INPUT THE DATA NEEDED BY MICROPEP * 

c * 


CRRA RARRRERKKKKRARARARARRRERREAKAEAAERARREREERARARERRRAEREKEER 


: , 
1010 call inputl(pc,pamb) 
write(*,40) 
40 format(1lx,'Are there any corrections of input for micropep?,y-yes 
*n-no ') | 
read(*,20)yes 
if(yes.eq.'Y'.or.yes.eq.'y') go to 1010 


CRARRARRRRKRKRERRRRRKKRAKKERERREKREAERRKARAERKERAARERERRERARERERERKRE 
Cc *® 


c¢ COMPUTE THERMAL EQUILIBRIUM, PROPELLANT CHARACTERISTIC, * 
C AND ETC. (USING MICROPEP, UNTILL GOOD RESULTS ARE REACHED) * 
Cc * 
CHR RARKRREKAKAKAEEKREKEKERRKEKKKKARREREKAREKAKKKKRKKRKRKRKKK 
S 
call micropep . 
write(*,50) , 
50 format(1lx,'Are the results good? y-yes neno ‘') 
read(*,20) yes 
if(yes.eq.'y'.or.yes.eq.'Y') then 
go to 1020 
else 
55 write(*,60) 
60 format(1x,'(1) Change the Basic Inputs of motor.'). 
write(*,70) 
70 format(lx,'(2) Change the Input data for Micropep.') 
write(*,80) 
80 format(l1x,' type 1, or 2, if 1 is chosen it is necessary') 
write(*,90) . 
90 format (1x, ' that Micropep shall be run once again!!!!') 
redd(*,100)i tg 
100 format(il) . 
if(i.ne.l.and.i.ne.2) then’ + 
write(*,*)! Error input!! 1, or 2 must be typed’ 
- go to 55 .- 
else 
if(i.eq.1) go to 1000 
if(l.eq.2) go to 1010 


endif 
endif 
RARER REREKERKERRERRRRAREREREERKREKRARAKRRREKRKS 
c * 
c CALCULATE THE LOSS PERFORMANCE & COMPUTE THE DATA FOR * 
c GRAIN DESIGN * 
Cc ¥ 


CRRA AAERREKRREKRERERERRRERAKE 
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C 
1020 call pfe(favg,tb,1,d,pc,cfc,ath,etacf,m,vp,vl,dpe,aport, 
*pfac,isp,cfth,webrqd,wf,tc,rrqd,ae,etac, xXx) 
write(*,50) 
read(*,20)yes 
if(yes.eq.'y'.or.yes.eq.'Y') then 
go to 1030 . 
else ® 
75 write(*,60) 
write(*,70) 
write(*,71) 
71 format(1x,'(3) Change nozzle data.') 
write(*,72) : 
72 format (1x, typerl, 2, or 3 ) 
read(*,100)i ; 
weiicit.leand.i.gt.3) then 
write(*,*)'Error input!! 1, 2 or 3 must be typed!‘ 
go to 75 
else 
if(i.eq.1) go to 1000 
if(i.eq.2) go to 1010 
if(i.eq.3) go to 1020 
endif 


endif 
CRRKRKKKKKKKRKKKKRKKKKKRKRKRKKKRKKKRKRKKKKKKKKKKKRKKKKKKKKKKKKKKKRKKK KS 


CHOOSE THE APPROPRIATE GRAIN DESIGN 

NOTE: OPTION 2, 3, 4, 5, 6, 7 ARE NOT AVAILABLE. INPUT 
YOUR OWN DESIGN BY SUBSTITUTE THE WRITE AND GOTO STATE- 
MENTS IN THE CORRESPONDING BLOCK 


qaqqagaaaa 
* t+ 4 + & 


CQRKKKKRKRKKRKAKRKRKKRKKKKRKRKKKKKKRKKRKKRKRKKKKKKKRKKRKKKKKRKKRKRKRKRKRKKKKKRK KKK 
1030 write(*,110) 
110 format(1x, 'Choose one of the following grain configurations') 
write(*,120) 
120 format(1x,'(1) Star grain’) 
write(*,130) 
130 format(1x,'(2) Optional grain ') 
write(*,140) 
140 format(1x,'(3) Optional grain ') ee 
write(*,150) 
150 format(1x,'(4) Optional grain ') a 
write(*,160) ape 
160 format(1lx,'(5) Optional grain ') ~* 
write(*,170).. < : 
170 format(1x,'(6) Optional grain !) 
write(*,180) 
180 format(1x,'(7) Optional grain ') 
write(*,190) 
190 format(ix, 'Type 1, 2, 3, «e---etc.') 
read(*,*)i 
pE(i.1t.lwor.i.gt.7) then 
write(*,200) 
200 forum tx (Error input! ! mwaececvpe 1,92, Zyuweoeo,7') 
go to 1030 
else 
if(i.eq.1) then 
rr5=d/2.0 
write(*,*)'Run "DESIGN MODE" ! 
1050 Call input2(pfac,aport,webraqd,rr5,1,tb,xx) 


vel 
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280 


290 


300 


310 


320 


330 
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call grains 
write(*,*)'Are the results are acceptable? y=-acceptable n-not 


*acceptable type y orn..' 


read(*,20)yes 
if(yes.eq.'n'.or.yes.eq.'N') then 
go to 1050 
else 
go to }055 
endif 
rr5=d/2.0 
write(*,*)'Run the "BURNBACK MODE" to get web vs. burning area’ 
call input2(pfac,aport,webrqd,rr5,1,tb, xx) 
call grains 
write(*,*)'Are the results are acceptable? y-acceptable n-not 


*acceptable type y oron..' 


read(*,20) yes 
if(yes.eq.'n'.or.yes.eq.'N') then: 
write(*,*)'Do you want to sry another star grain? y/n : 
read(*,20)yes 
if(yes.eq.'n'.or.yes.eq.'N') then 
goto 1030 
else 
go to 1050 
endif 
else 
go to 1060 
endif 


else 


if(i.eq.2) then 
write(*,280) 
format(1x, ‘Input your own grain design package’ ) 
go to 1030 
else 
if(i.eq.3) then 
write(*,290) 
format(1x, 'Input your own grain design package') 
go to 1030 
else 
if(i.eq.4) then 
write(*,300) 
format(1x,'Input your own grain design package') 
goto 1030 . 
else o 
- df(i.eq.5) then 5 
write(*,310) 
format (1x, "input yours own ‘grain design package') 
go to 1030 
else : : 
if(i.eq.6) then 
write (*,320) 
format(1x, ‘Input your own grain design package') 
go to 1030 
else 
if(l.eq.7) then 
write (*,330) 
format(1x, ‘Input your own grain design package') 
go to 1030 ee 
endif 
endif 
endif 


U2 


me 


m= «tee ge ee oe 
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endif 
endif 
endif 
endif " 
endif 
e 
CRRA KEKEREKEREKEKEEEKKERREKKEKEEKERKEEREKEEKERKEEER 
¢S ® * 
c INPUT THE DATA FOR THRUST & PRESSURE VERSUS TIME PROFILE* 
c¢ CALCULATION * 
c * 


CRHRAKKKKREAKERKKKEEEKEEEEAKEKEKEEKEKKKEEKKEKKKEKKKEKKRKRKKRKKEKKES 


1060 write(*,*)'The following processes is to calculate a "THRUST VS. 
* TIME" and “PRESSURE VS. TIME" profiles’ 
va=(d/2.0) **2.0*3.14159*1 
rrqd=webrqd/tb 
write(*,106) rrqd 
106 format(1x, ‘Required burning rate for new web, IN/SEC.. ',f£6.4) 
1070 call rektin(ath,ae,va,tamb,pamb,pc,etacf, rrqd, pik, etac) 
write(*,*) 'Any change EOL the input? y-yes, n-no 
read(*,20) yes 
if(yes.eg.'y'.or.yes.eq.'Y') then - 


go to 1070 
else 
go to 1080 
endif 
Cc 
QRRRRARARRAKEKKEKKEKREKEKKERERKRKERARKRKRAEKRRKAKKKKKRKRKKRKKKRKRKRKKKRKKK 
e * 
c¢ COMPUTE THE PRESSURE VS. TIME & THRUST VS.TIME PROFILE * 
Cc OF ROCKET MOTOR * 
Cc se 


CRRA EKEKEKKKKKKKRKEKRKKKKKEE 
c 
1080 call rocket 
¢ 
write(*,340) 
340 format(lx, ‘Are the results acceptable? yryse, n-no......') 
mege ("720) yes 
if(yes.eq.'y'.or.yes.eq. 'Y') then 
go to 9999 
endif <i 
write(*,350) = 
350 format(1x,'Do you want to proceed?” YrYES, N~“NOw ce cevcccee’) 
read(*,20)yes” .- ~ 
if(yes.eq.'N'.or.yes.eq.'n') then 
write(*,*)'Terminated on the request of user ' 
go to 9999 
endif 
1040 write(*,360) 
360 format(1ix, 'Suggest to change the following data') 
write (*,380) 
380 format(1x,'(1) Change the Grain Design. ') 
write(*,390) 
390 format(lx,'(2) Change the Ingredients of the propellant. ') 
write(*,400) 
400 format(1x,'(3) Change the Basic Design Requirements. ') 
write(*,*)'(4) Re-run the "ROCKET"! 
write(*,*)' Type 1, 2 or 3 for selection' 
read(*,100)i 
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if(i.me.lcand.i.ne.2.and.i.ne.3) ‘then 
write(*,%*)' Error input you have to type l, 2 or 3! 
go to 1040 
else 
if(i.eq.1) then 
go to 1030 
elise 
if(i.eq.2) then 
go to 1010 
else 
if(i.eq.3) then 
go to 1000 
else 
if(i.eq.4) go to 1070 
endif : 
endif 
-ndiz 
endif 
CARRRRAAAARARRERERRARRERRRRRRRERRERERRRRRERRRRRRAKRRRKRKRRRKKKRRRAE 
Cc ® 
c INPUT THE TRAJECTORY DATA (FUTURE WORK) * 
c 2 
CRRKRKRRRRKRKRRRRRKRKRRRKRKKRKRKRRRKRRRKRRKKRRRRKRRKRRRRRKRRKRRRKRRRRRKRKRK 
c 


e call flyin 

Cc 
CRRRERERERRRRRRRRERRRRERRRRERRERRRERRRRRRRERRERERREKRRRAKRRRRRRER 
c * 
c¢ CALCULATE THE TRAJECTORY OF ROCKET MOTOR (FUTURE WORK) * 
G * 


CRRREREERERERARRARARERERRARARERRARERRRRRERRARRRRRERRERRRARRRKRRRRRKRRKER 


c 
c call flyit 
9999 end 
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subroutine Input0(favg,tb,1,d,meop,pc,r,pamb,tamb,h,fmin,tpmax, 


*tpnom, tpmin, pik) 
QURRKRKRRKRKRKRKKKRRKRKRKKKRKRKRKKKKKR KR KK KKK KKK KKK KKK RRA RRR 
c * 
¢ THIS SUBROUTINE IS FUNCTION AS A BASIC INPUT FOR SOLID * 
¢ PROPELLANT ROCKET MOTOR DESIGN * 
Cc 6 * 
‘ohEeEE EER ERLE LE LEE LEE EE EEE PELE RR E EEE EET ELE EERE LER EE EEE SE ES 
e 
c variable definition: 

G eukkkkkaee (ballistic performance) *****xkees 
Cc favg >; average thrust (F12.6) 
c tb : burning time of motor (F12.6) 
e fmin ¢: minimum thrust (F12.6) ’ 
c kekkhkhkkak (mission/vehicle) RkekkaRRKKRKRKR RRR AS 
S 1 : motor length (F12.6) 
Cc a : motor diameter (F12.6) 
Cc meop ; maximum expeted operating pressure (F12.6) 
Cc va : motor volume avaiable for propellant (F12.6) 
S tpmax, tpnom, tpmin 
e pmaxtp : meop * (l-reproducibility tolerance) (F12.6) 
e pe : norminal chamber pressure (F12.6) 
S rptol ; reproducibility tolerance (F12.6) 
e r : burning rate (F12.6) hii 
Cc web : web of grain (F12.6) 
c pik : temperature sensitivity of propellant(F12.6) 
e h : altitude(F12.6) 
real favg,tb,fmin,1,d,meop,va,dltp,pmaxtp,pc,rptol,r,web,h 
c open a file which can hold the basic input data 


open(2,file='input0.dat',access='sequential',status='unknown') 

e input basic data 
write(*,*)'Al1l input data in decimal format’ 
write(*,*)'Specify the average thrust, LBF........ ' 
read (*,10) favg 

10 format (f12.6) 
write(2,10) favg 
write(*,*)'Input burning time of motor, SEC....... ' 
read(*,10)tb 
write(2,10)tb 
write(*,*)'Input minimum thrust required, LB...... ' 
read(*,10) fmin ae 
wate (*;m) ‘Imput motor length, INCHES. ....ccsvisiocee * 
read(*,10)1l a 
write(2,10)1 
write(*,*)'Input motor diameter, INCHES. .... 0.0005 
read(*,10)d . 2 
write(2,10)d ? 
write(*, *) ‘Input maximum expected: ' 
Wriuee(*,*)'aperating pr@ssure, ESI. . ccc cece vcese cum ! 
read(*,10)meop "4 
write (2, 10)meop 
Wren ey)” TNDUC Tomax, Fecccccescccecccccccccccce ' 

‘ read(*,10) tpmax 

write (2,10) tpmax 
Wire (*,*) Ineat TorGm, ders «os ccs 006 6 0 elcicc siellviles 00 
read(*,10) tpnom 
write(2,10) tpnom 
Lot * we rmOut Tm, .Puiuiuensie 0 «leis « oisie sMeieleKeueee ee 
read(*,10)tpmin 
write(2,10)tpmin 
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write(2,10)dltp 

write(*,*)'Input reproducibility tolerance, cian : 
read(*,10) rptol 

write(2,10)rptol 

rptol=rptol/100.0 

pmaxtp=meop*(1.0-rptol) 

rptol=rptol*100.0 

r=(d/4)/tb 

write(*,20)r 

write(2,10)r 

format(1x, ‘Estimated repuired burn rate.......ecee- ',£12.3, 


*' INCH/SEC') 


write(*,30) 

The following calculation of ambient ommaeeiee and préssure at. 
certain altitude is directly copied from NWC Technical Memorandum 
4757 subroutine atmstd based on U.S Standard Atmosphere, .1976. 
format(1x,'Input the required design altitude of rocket motor in F 


eT eeeoveeveeveeveeeeeeeeeeeeeeeeeeeereeoeeveenenene 


xt 


read(*,10)h 

write(2,10)h 

Lif (h-36152.)1,1,2 

Camb=518.67-0.003559969*h 

ro=0.00237696*exp(-(0. oaeE EQ 2a Gane 001)+0.000132317* (h*.001) **2) ) 
go to 7 

if(h~-65825.)3,3,4 

tamb=389.97 
ro=0.00237696*0.29697*exp (-0.04771923* (h-36152.)*.001) 

go to 7 

if(h=105511.)555,6 

tamb=354.142+0. 0005442926*h 
ro=0.00237696*0.0718594*exp(-(0. Oe te )*.001 


&=0.000034655*((h-65825.)*. anes 


go to 7 

tCamb=251.4114+0.00151794*h 

ro=0.00237696%*0.0107993 *exp (= (h-105511. )}/21320.) 
pamb=11.9203*ro*tamb 

write (2,10) pamb 

write(2,10)tamb 

end of calculation of tempressure and Genp eneaire at deaean alt. 


write(*, *) ‘Based on available requirements, make an indie 
write(*,*)'! selection of propellant. Then-input the temperature! 
write(*,*)'sensitivity of the propellant, %7r... eee 


read({*,40) pik : 
write(2,40)pik 3 
format (f12.6) ¥ 


pik=pik/100.0 ; 
pc=pmaxtp*exp (pik* (tpnom=-tpmax) ) 
pik=pik*100.0 

write(2,10)pc 

close(2,status='Keep') 

return 

end 
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subroutine inputl(pc,pamb) 
CRARRRRARKKRRRERRRKRKKKRRRAKKERKRKRRKRARRERRRRRRARKKRKRRRRRKRKKEE 


open(3,file='input.dat!,access='sequential',status='unknown') 
c input the case name from the keyboard 
write(*,*)' Start to input the selected propellant data a 
| write(*,*)' ! 
! write(*,1) 
| a format (1x, 'Specity the case name (format Al0O)...... '") 
i 


c * 
c¢ THE PRIMARY PURPOSE OF THIS SUBROUTINE IS TO ESTABLISH * 
c¢ The DATA FILE WHICH CAN BE ACCESSED BY MICROPEP SUB= * 
c TINE * 
e * 
CRERRKREKEREKRERRERARERARERREREKKERKRRKKKREKRKKKRRKKKRKKRKKKEKRKKKKRKE 
ec 
i common /micrp/ a(12,12),kr(20),amat(10,12) ,jat(12),in,is, fie(10,6) 
1,ie(10,6),alp(12) ,w27,n,dah(10) ,rho(10) ,wate(10) ,wl(6),w43,1ig,np, 
: 2vnt(201), w47, name,ser, floor, itag(100) ;wing (10) 
? common/chara/block, je,aspec,specie . 
! character*8 specie(200) a: 
character*30 block(10),ofile 
character*2 je(10,6),aspec(12) 
common/moon/tstest,te,irun, iounit, iend 
| character*30 blok(10) 
character*10 user ~ 
dimension jie(10,6) 
integer nic,niu,irun,alt,kr 
aa denexp 
es - 
C open the file as formatted sequential as unit 3 


read(*,'(a)')user - 
write(3,10)user 
10 format.(1al10) 
irun=2 
c¢ input the number of ingredients, supplied ingredients, and runs 
write(*,*)'Input the number of ingredients? (format I5)....' 
read(*,*)nic 
write(*,*)'Input the number of user suppplied ingredients (I5)...! 
read(*,*)niu 
WEice(a,25)nic,niu, irun 
20 format (2i5) 
25 format (3i5) 
c¢ Input density exponent which is used in conjunction with control 
© Optien 3. 
write {*,*)'Input the density exponent: '! 
write(*,%*)'For tactical missiles or first stages use 1.0' 
write(*,*)'Use 0.7 for second stages and 0.2 for third stages' 
Vatmmue, ©) ' (format FLO.O)...ssseccsvecesctoue 
read(*,40)denexp 
write (3,40)denexp 
0 format (f10.5) 


= — = 


input control option store it in kr(20) 
print out the options on the screen 
write(*, fa)! )'There are five options listed below for your' 
write(*,*)'selection. Type l, 2, 3, 4, 5---' 
ee ae a ouren input for design is 3! 
write(*,'(a)')! 1. For chamber and exit calculations. ' 
f 
e 


qa0a a & 


write(*,'(a)')' 2. For chamber only calulations.' 
write(*,'(a)')' 3. For boost performance and nozzle design.' 


ey 
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write(* ita). ) ‘ 4. For equilibrium calculations at specified T' 
write(*,'(a)')' and P. Use T in Kelvin at Pe location and P' 
write(*,'(a)')! in psia at the normal Pc location.' 
write(*,'(a)')'! 5. Other options ( Check README for special' 
write(*,'(a)')'! alternatives).' 


write(*,*)'Use format Il ' 
read(*,4S5)alt 
45 format (il) 

if (alte. ltwS) then 
kr(2) = 0 
kr (4) 0 
kr(5) 0 
kr (6) 0 
kr(8) 0 ; 
if (ait.eq.1) then 


en 
kr(1) fe) 
kr (3) 0 
kz 3 0 
else 
if (alt.eq.2) then 
kr(1l) = 1 
kr(3) = 0 
kr(7) = 0 
else 
if (alt.eq.3) then 
kr(1) = 0 
kr(3) = 1 ° 
kr(7) = 0 F 
else ae 
kr(1) = 0 
kr(3) = 0 
kr 7). 2: 
endif 
endif 
endif 
write(3,50) (kr(i),i 
50 format (811) 
else 
write(*,'(a)')' Input your requirements for calculation--' 
read(*,60) (kr(i),i = 1,15) 
write(3,60) (kr(i),i = 1,15) 


1,8) 


60 format (1511) 

endif 
c input the serial number of ingredients found in JANNAF. DAT table 
Cc 


write(*,70) nic 2 : 
70 format(1x, 'Input : ,i2,' ‘ingredients needed for calculation") 
do 100 iwi,nic . oo 
write(*, 80) 4 


80 format(1x, 'Serial # of number ',1i5,' ingredient(I5)... ) 
read(*,90) itag(i) 
90 format (i5) 


100 continue 
write (3,120) (itag(i), i = 1,nic) 
120 format(10i5) 
c Input user supplied ingredients 
if (niu.gt.0) then 
do 170 i = l,niu 
write(*, ' (a)! )'Specify ingredient name use caps lock (A30) 
read(*,130)blok(i) 
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130 
140 
150 
160 
170 


172 
175 
180 
190 


210 


220 
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write(*,'(a)')*Specify composition of ingredient max = 6! 
write(*,'(a)')' use (6(1I3,A2))' 
write(*,'(a)')! PRARAPHHAAHRP AREA HARE EE RRA EE RH! 
read (*,140) (jie(i,1) ,je(i,1) ,1=1,6) 
write(*,'(a)")? Specify heat of formation, and density (f£5.0,f6 
*.0) 
read(*,150)dh(i) ,rho(i) 
write(3, 160) blok(i), (jie(i,1),je(i,1),1=1,6),ah(i),rho(i). 
format(a30) 
format(6(i3,a2)) 
format(f£5.0,f£6.0) 
format (a30,6(i13,a2) ,£5.0,£6.0) 
continue 
else 
write (*,'(a)')' No user supplied ingredients. ' 
endif 
if, (igun. 1lt12)ethen 
irun = 2 
endif 
status=0 
do. 220 i = 1,irun 
if(status.eq.0) then 
w1(5)=250 
-" write (*,'(a)')' sneer the desired chamber pressure (F1l2.4) ' 
read (*,180)w1(5) 
wl (6)=pamb 
_ write (*,'(a)')' Specify the exhaust pressure (F12.4) ' 
read (*,190)wl1(6) 
ni=nic+niu 
do 175 j =1,ni 
write(*, 172) itag(4) 
format(1x,'Input the weight % of =nameaaent' . iS, "Gero. 8). . 1") 
read(*,190)wing(j) 
continue 
format(f10.4) 
format(£10.4) 
write(3,210)wl1(5),w1(6),(wing(j),3=1, oe) 
format (12£10.4) 
status=l 
else 
w1(5)=pe 
write(3,210)w1(5),wl(6), (wing(3), =i, 10) 
endif: a 
continue ; 
close (3) 
return 
end 


uo 


APPENDIX 


subroutine pfe(favg,tb,1,d,pc,cfc,ath,etacf,m,vp,vl,dpe,aport, 
+pfac,isp,cfth,webrqd,wf,tc,r,ae,etac, xx) 
CRRAAARAARAARARRARARRRRRARRRRARRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRR 


Cc ® 
c¢ THE MAIN PURPOSE OF THIS SUBROUTINE IS TO CALCULATE THE #* 
Cc PERFORMANCE LOSS AND THE DATA FOR GRAIN DESIGN * 
Cc . x 


CRARAARARAARRAAARAARAARRERRRRRRRRRARRRRRRARRRRRRRRRRRRRRRRRRRRRRE 
C variable definitions: 


c favg: Average thrust from basic input 

e cLen: thrust coefficient from micropep 
Cc cre -; corrected cf 

c ewe: chamber pressure from basic input 
c qd : motor diameter from basic input 6 
c i ; motor length from basic input 

C athth: theoretical throat area 

Se ath: throat area 

¢ atht ; tempory result of throat 

c athe: theoretical throat diameter 

c etacf: thrust coefficient efficiency 

S Gia See for convergence test of throat size 
G va : vol avaiable for propellant 

e ED : burn time 

c vp : propellant vol. 

S vl : volumetric loading factor 

e aport: port area 

c dpe ; equivalent circular port diameter 
c pwet ; wetted per_meter 

C 5 ; throat-to-port area ratio 

c rHop: propellant mass density 

C mw : molecular weight of propellant 

c mf : mole fraction of condensed phase 
S webrqd: required web 

c wf a web fraction 

c prac perimeter factor 

e character*8s0 fname 

C 


real favg,cfth,cfc,pc,d,1,athth,ath,dthth, etacty ait, tb, Vp, Vleearae 
+aport,dpe, rhop7a,nn, rstd ,rpc,web »webrad, wf,m, pwet, pac, isp, ispm, at 
+ht,mw,t¢ aspen, 1a, mf 
Character*l yes 
c The purpose of this common instruction is to transfer aai(oropelaeee 
c density) from MICROPEP: SUBROUTINE PUTIN(LE) 
common /micrp/ a(12,12),kr(20) ,amat(10,12),jat(12),in,is, 
lfie(10,6),1e(10,6) ,alp(12) ,w27,n,dh(10) ,rho(10), 
2wate(10) ,wl(6) ,w43,ig,np,vnt(201},w47,name,ser 
3,floor, itag(100) ,wing(10) : 
C The purpose of this common instruction is to transfer the CHAMBER PRE- 
C SSURE & MOLECULAR WEIGHT from MICROPEP SUBROUTINE OUT. 
common /tcf/ tfc,rmwtc 
c Transfer SPECIFIC IMPULSE, CHARACTERISTISC VEL., CODFF. OF THRUST & 
C OPTIMUN EXPANSION RATIO from MICROPEP SUBROUTIN DESING. 
common /socf/ spi(2),cf(2),c0ex(2) 
c calculate the theoretical throat area 
, pi=3.14159 
c transfer cfth from "MICROPEP" 
cfth=cf (2) "aoa8 
Cc initiabizercc. 
cfc=cfth 
athth=favg/ (cfth*pc) 


80 


‘ 
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c initialize ath=z=athth 


ath=athth 


co calculate theoretical throat diameter 
dthth=saqrt (ath/pi) *2 
write(*,888)cfth,athth,dthth 


888 


fermac ix, Cren™= ,f5.3,' 


c calculate the etacf 
write(*,*)'Input the nozzle half angle! 


write(*,*)'alpha, 


read(*,*)alpha 


- pi=3.14159 


alpha=(alpha/l 


80.0) *pi 


Athth=',f6.3,' Dthth=', 


write(*,*)'Input the nozzle exit half angle’ 
Wreee @ ) Chetaex, DEGREES (F5.2) .ccccccccvceves 


read(*,*)theta 


ex 


thetaex=(thetaex/180.0) *pi 


write(*, *) ‘input radial erosion rate of the throat in IN/SEC! 
write(*,*)'assumed to start at time zero (FORMAT 0.XXXX)....! 


read(*,*)erosr 


write(*,*)'Input the mole fraction of condensed! 
write(*,*) ‘phase in (MOLE/100 JG) eee cccvavccccvececs 


read(*,*)mf 


write(*,*)'Input the material of nozzles! 


write(*,*)'l-steel nozzle, 


read (*,*)nozl 
atht=ath 
do 9 i=1,50 


e 


C calculate the-loss mechanism 
call loscf(etacf;tb,ath,pc,alpha,thetaex,erosr,mf,nozle) 
eGvecorrect the trust B26 ee 


cfc=cfth*etacf 


c calculate the corrected rene area 


8888 


ath=favg/ (cfc* 


pc) 


write(*,8888)etacf,cfc,ath 
3,° Cher’ #fa2.3," . Atha’ ,f12.3) 
c¢ calculate the convergence athth and ath 
dif=abs ((atht-ath) /atht) *100.0 


format (1x, ‘Eta 


fe (aetegt. 1.0 
atht=ath 
else 
go to 91 
endif 
continue 
va=pi*(d/2) **2 
etac=0.93 
vol=pi*(d/4) ** 


efe™,f12. 


) then 


*) 


2*1 


O-ordinary nozzZle..... 


“=r 


c Density and molecular weight com from UMTCROreP" 


98 


99 


87 


rhop=w43 
mw=rmwtc 
tc=tfe 
ispth=spi (2) 


b 


Gere 6a\rhonam. tC, 1Speh 


format(' Rhop=', 

a) 
isp=ispth*etac 
m=favg/isp 


bp aeele2 3,' 


f*etac 


vp=favg*tb/ (isp*rhop) 


l=vp/va 


write(*,87)VP,VA,VL 


format(' VP= 


price, ° 


Mw=',f12.3,' 


VA=',F12.1,' 


8 1 


£5.3) 


DPCEME SME 2 \xccausevesccederes ' 


Tc='! peude-1, : 


VL=!,F12.3) 


Ispth=', 


f12.2 


-— 9 eee ee 
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dpe=sqrt(1.0-vl) *d 
aport=(pi/4) *d**2*(l-vl) 
tres=(aport*1) *pcexmw/ (m*#12*1545% (tc+459.6) ) 
tres=tres*1000.0 
write(*,88)tres 
format(' RESIDENCE TIME IN BS 922 "je Loe) 
write(*,89)etac 
format(’ COMBUSTION EPFICIENCGY (etac).-...ccceee. place or 
write(*,*) 'Change etac? y-ye@S, NenO.ecccccccscccce ! 
read(*,20)yes 
if(yes.ge.'y'.or.yes.eq.'Y') then 
write(*,*) "Inpit GtAC. . oso oo s-cvsjsisisisuesersueier terest terete 
read(*,30)etac : 
vol=aport*] : 
go to 99 
endif 


c calculate the noceate exit area 


ae=ath*oex (2) 


c¢ calculate the required perimeter factor 


2. 


22 


25 


24 


ig 


18 


write(*,19)pc | 
format(1x, ‘Input the: c* with Pecnome"/)%10.1,' (F10.4)..... ") 
read(*,*)cstar 
XX=(pc*ath*32.17)/ (cstar*rhnop*1*3 .14159 *dpe) 
write(*, eens the estimated required web ene (F10. 4)! 
write(*,*) 'suggest initial guess D/4......... 
read(*,*)webrqd : 
r=webrqd/tb 
pfac=xx/r 
write(*,2l)r 
format(lx, 'The required burn rate =',f6.4) 
write(*,22)pfac 
format(1x,'The required perimeter factor=',f6.3) 
write(*,25)xx me 
£LOormat (1x, 'pfage' m25.4,'/r') 
write (*,24)tb 
format(1x,'tb =',£7.2) 
write(*,*)'Are these acceptable? y-yes ne-no....... ! 
read(*,20) yes 
if(yes.ge.'y'.or.yes.eq.'Y¥') then 

pwet=pfac*pi*dpe . 

wf=webrqd/ (d/2) 

j=ath/aport 

ld=1/d - re 

go to 23 a 
else 

go to 18 
endif 
format(al) ; , 
format (f12.6) 
write(*,40)cfc 
format(lx, 'THRUST COEFFICIENT aiiicic << css ss ccueeee. | eee) 
write(*,50)etacf — 
format (1x, "ETACF. . «ss oc so cieieleleie cs 6 clcleleieeetets ctctere «(enum 
write(*,60)isp 
format(1x,'SPECIFIC IMPULSE, LBF SEC/LBM.........',f12.3) 
write(*,70)m 
format(1x, 'MASS FLOW RATE, LBM/ SEC « sieusieis sousieie 's =) peeeerre)) 
write(*,80)vp 
format(1lx, "PROPELLANT VOLUME, CU IN.ccccccccsccse 'p,LlZ-l1) 
write(*,100)ath 


wr ws 
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fommut(ix, "NOZZLE THROAT AREA, SQ IN.....-.seeceee', £1223) 
write(*,105)ae 

format(1x, "NOZZLE EXIT AREA, SQ IN.....e.eeeeeee', £123) 
write (*,110)dpe 
format(1x, ‘EQUIVALENT PORT DIAMETER, INe.cccccceee', £12.33) 
write(*,120)aport ; 
format(lx, 'PORT AREA, SQ IN a a 2 Fo i Ph 
write(*,130)pfac 

format (1x, "PERIMETER FACTOR... .ccccccccesccccccccce', £12.3) 
write(*,140)webrqd 

Poimieeeloe Re OULD WED ET Lileoclsiecccecceccc ccc cn cece stlec3) 
write (*,150) wf 
Commack, ROOULREDAWED BPRRCTION. ...cccerccccccvces yy fiee3d) 
write(*,90)vl 

format(1x, 'VOLUMETRIC LOADING... cccccccccccccseee',f12.3) 
write(*,190)la . 
format(lx, 'LENGTH TO DIAMETER Drammen. .... 00. et Pike. 3) 
write(*,170)tres 

BoLmemct@eer res teres LIME, MSEGasiwc cc cc cc ccc cc ccce 'p kbar) 
write(*,180) 4 

format(lx,'THROAT TO PORT AREA RATIO (J FACTOR)...',£12.3) 
write(*,200)r 

format(1x, "REQUIRED BURNING RATE, IN/SEC.......! RS). 3) 
return 

end 


te 
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subroutine loscf(etacf,tb,ath,pc,alpha,thetaex,erosr, 
*mf,nozle) 


CRRRKRRKRKRKRRRRRRRRRRRRKRRKRKRKRRRRRKRRKRKRKRKRRKRKRRRKRKRKRKRRRKRKRRRKRKKRKRK RK 


aAnAaAnAA|A 


*® 
THIS SUBROUTINE CALCULATE THE NOZZLE LOSS COEFFICIENT * 
(THRUST COEFFICIENT) USING THE METHODS PROPOSED BY AGARD * 
PROPULSION AND ENERGETICS PANEL FOR METALLIZED PRO= * 
PELLANTS (Reference: Performance of Rocket with Metal=- * 
lized propellants, AGARD-AR-230) a 

oe 

& 


CRARRRRKARKRKRRRKRKKRRKRRKRRKRRKRRRRKRKRKRKRKRKRRKRRKRKKRRRKRKRRKKKRRRKR RRR RRR 


qaaa 


0 Q 


qagqagaagagagagaagagaagaagnagnaagnaaangagqaagqaaaaa 


Defination of varibles'in this subroutine 


etacf: nozzle loss coefficient 
divs: percent divergence loss ’ 
kin : percent kinetics loss 
bl : percent boundary layer loss 
tp $ percent two-phase loss. 
sub 3 persent submergence loss 
eros : percent nozzle erosion loss 
astar: nozzle entrance/ nozzle throat area 
Ss : submergence length/ length of internal motor 
mf : mole fraction of condensed Jee in mole/100grm 
C1,¢C2Z,7¢63,¢%,Co - 
constant in emperical relation 
dp : diameter of particle 
pe $ chamber pressure 
ispf : Isp of frozen equilibrium 
isps : Isp of shifting equilibrium 
nozle: steel nozzle = 1, ordinary nozzle = 0 
tb : burning time of motor 
at : throat diameter 
expn : nozzle expansion ratio 
ispi : Isp with ideal (at design point) nozzle expansion ratio 
ispm : Isp with mean nozzle expansion ratio 


integer expnu,expnl 
real etacf,div,kin,bl,tp,sub,eros,astar,s,mf,dp,pc,isps,ispn, os 
+dt,ispu, ispl, erosr 
common /scratc/ plot(5,100) 
common /socf/ spi(2),cf(2),0ex(2) 
calculate the % of nozzle divergence loss 
div = 50 * ( 1 - cos ((alpha + thetaex) / 2 )) 
ispf is the frozen and isps is the shifting specific impulse from 
micropep run ‘at nominal pc 
Transfer specific impulse data from MICROPEP.° 
ispf=spi(1) ake. ? 
isps=spi (2) + : 
Calculate-the kinetic loss 
kin = 33.3 * (1 - ispf/isps) 
calculate % boundary layer loss 4 
Steel nozzle (1), Ordinary nozzle (0). Enter 0 or 
if (nozle.eq.1) then 
cl = 0.00506 
c2 = 0.0 
else . 
cl = 0.00365 a 
c2 = 0. 000937 ; 
endif 
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expansion ratio of the nozzle for shifting equilibrium from 
micropep with exit pressure=ambient pressure at design altitude 
read (*,70)expn 
expn=oex (2) 
pi=3.14159- 
dt=sqrt(ath/pi) *2 
bl =(cl*pc**0.8/ (dt#0. 2)) *(1+2*exp (-c2*pc**0. B8*tb/ (dt*0.2)))*(1+0. 
$016*(expn-9.0) ) 
calculate the two-phase flow loss 
the mole fraction of- condensed phase is come from "MICROPEP" 
mean particle diameter calculation by using empirical relation — 
dp = 3.39*dt*#*0.4692 , 
if (mf.1t.0.09) then . ; 
mf = 0.09 : 
endif 
ti(ct.1t.1.0) then 
c3=9.0 
c5=1.0 
c6=1.0 
else 
if(dt.le.2) then 
| c3=9.0: 
c5=1.0 
c6=0.8 as ; 
else 
ncap. lt. 4) then 
c3=13.4 
c5=0.8 
c6=0.8 
else 
if(dp.le.8) then 
c3=10.2 
c5=0.8 
c6=0.4 
else 
c3=7.58 
c5=0.8 
c6=0.33 
endif 
endif 
endif 
endif oo 
c4=0.5 = 
tp=c3*mf**c4*dp**c5/ (pcox#0.15*expn**0.08*dt**c6) 
calculate the submergence loss(NOT USED IN PRELIMINARY WORK !!) 
the nozzle is assumed to be of external design with no submergence 
loss. s=length of submergence / length’ ofinternal motor. 
astar=(pi*d**2/4) /ath 
sub=0.0684*(pc*expn/astar) **0.8%*s*#0.4/dt**0.2 
default to 0.0 
sub=0.0 . 
calculate the nozzle erosion loss 
calculate the ispm (specific impulse at mean expansion ratio) 
df=dt+erosr*tb*2 
dbar=(dt+df)/2 
expnm=(dt**2/dbar**2) *expn 
expnl=int (expnm) 
Transfer expansion ratio of nozzle and specific impulse from MICROPEP. 
ispl=plot(4,expnl) 
expnu=int(expnm) +1° 
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ispu=plot(4,expnu) 
ispm=(expnm-expnl) / (expnu-expnl) *(ispu-ispl)+ispl 


Cc erosion loss calculation 


100 
110 
120 
130 
140 


150 


160 


eros=(l-ispm/isps) *100 

write(*,100)div 

format(1x, 'PERCENT DIVERGENCE tBSS%.... Mes ee ',£6.3) 
write(*,110)tp 

format(1x,'PERCENT TWO PHASE FLOW LOSS........++++2-- ',£6.3) 
write(*,120)bl : 
format(1x,'PERCENT BOUNDARY LAYER LOSS. ..eeseeeeeeeee ', £6.3) 
write(*,130)kin 

format (1x, 'PERCENT KINETIC LOSS...-.cescccsccsccsccee .',£6.3) 
write(*,140) sub 

format (1x, ‘PERCENT SUBMERGENCE LOSS..cccccccccccssccs mmptoOes) 
write(*,150)eros 

format (1x, ‘PERCENT EROSION LOSS .ccccccccccccccccsccese ‘, £623) 
etacf=1-(div+tptbl+kin+sub+eros) /100.0 

write(*,160)etacf 

format (1x, "THRUST COEFPFECIEIT EFPEGREN Gia. «<scs o os cam, BOsa) 
return 

end 


ve 


ay Tod 
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subroutine input2(pfac,aport,web,rr5,lgrain,tb, xx) 


CRRRARKKKKKKRKRKRRKRKRKRKRKRKKRKRKRKRKRKRRKRKKRKKRKRKKRKRKRKRKRKRKRKRK KKK RRR 


Cc. 
Cc 
Cc 
Cc 


* 


THE PRIMARY PURPOSE OF THIS SUBROUTINE IS TO ES=- * 
TABLISH THE DATA FILE FOR THE CSD "GRAINS" PROGRAM * 
* 


CERRRRKRKRKRKRKRKRKRKRKRKRKRKKRRRRKRRKRKKRKRKRKRKRRKRRKRKRRKRKRRRKRKRRRKRRRKRRRRRRRRREK 


pL. 
14 


13 


character*80 tape5 
character*l yes 
real lgrain, 1lspoke,kspoke,RR5 
integer nmax,nspoke,ibb,iprnt 
ae ene following will be the data input for ' 
write(*,*)'spoked star grain designs with radial auch : 
open(4, file='grainin.dat',status='unknown') 
write(*,*)'Run DESIGN-1 mode or BURNBACK=0? Type 1 or fo... 
read(*,*)i 
if(i.eq.1)then 
write (*,40)web 
=web/tb 
pfac=xx/r 
write(*,13)r 
Seemat (ix, "Requiredeourningwrate,IN/SEC. .cccccccccccce ',tl2.3) 
write(*,50)pfac 
write(*,*)'Do you want to Change the web?..cy/Meeecsccccee ' 
read(*,5)yes 
if(yes.eq.'Y'.or.yes.eq. 'y')then 
write(*,*)'Input the required WDM NicicieMeNI sc cleletstcs ccee 
read(*,*)web 
go to 14 
endif : 
Were se) MAX NUMDer Of SPOKeS?s(I3)iccenescccscesscccce. ' 
read(*,*)nmax 
oece wma Initial number of SPOKES?» (13) o0..+-esccismeccevee * 
read(*,*)nspoke 
ibb=1 
iprnt=0 
write(4,10)nmax,nspoke,ibb,iprnt 
write(*,30)aport 
30 Beumde Cre RECQULECGMCOLE SECA, SO Dia. ..cewesescsscececee ',f12 
*.3) 
write(*,40)web 


_ 40 Pormae (lx, 'Interalwebwenicknes sw aTlIN. ..ccccssiesicccecesee ', £12 


*.3) er 
write (*,50)nfac 
50 format (1x, 'Required perimeter factor Soossccooo0ol d+ Cee oi 
*.3) 
write(*,*) 'Input Kspoke - Abstract multiplier on ! 
write(*, *) 'PHI2 - if KSPOKE=0.0 PHI2=0.0 where: 
write(*,*)! Ee Ne CS pee Ee Bee COs ae) (0.0-1.0)e... ' 
read(*,*)kspoke 
write(*,60)rr5 
60 format ( ise, "Gain CUEEE PadlUS) LNCHES. «2.6 cccce ss cscs ccces', £12 
*.3) 
write(4,20)aport,web,pfac,kspoke,rr5 
write(*,*)'Input initial grain inner radius, INCHES (F7.3) ' 
read(*,*)ril 
write(*,*)'Input initial value of SR1 - spoke! 
write(*,*)'top corner radius in IN (Default=0.25)........-. 
read(*,*)sril 
write(*,*)'Input initial value of SR2 - spoke’ 
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write(*,*)'base coner radius in IN (Default=0.25)......... ' 
read(*,*)sri2 
write(4,20)ril,sril,sri2 
write(*,*) 'Input increment on SR1l in IN = cannot be 0.0... ' 
read(*,*)dell 
write(*,*)'Input increment on RR1 in IN - negative value.. 
read(*,*)del2 , 
write(*,*)'Input increment on NSPOKE Default = 1.0........ ' 
read(*,*)del3 
write(*,*)'Input increment on SR2 in IN = cannot be 0.0... ! 
read(*,*)del4 
write(4, ee: del2,del13,del4 
write(*,*)'Input tolerance on grain angle sum total - rt 
write(*,*) "of 360. (0.0=170)...... co cscmmene + < cc ccs Sides c. <sen 
read(*,*)cirtol 
write(*,*)'Input tolerance on port area - % of Aport ' 
Wraice (#7 (0.01.0) «+000 000 000 60¥ ee 5 6 6g ele n citi iinet nreenenCneEE ' 
read(*,*)pttol 
write(*,*)'Input tolerance on perimeter - % of perimeter.. 
wEeite (* 9%) (0. 0=1.0) uence «0 000s © silewoneuelblie Wels o ciceet tetas nnn nmmaC a: 
read(*,*)pertol 
write(4,20)cirtol,pttol,pertol 
write(*,*) 'Input Max Of SR1 LWP CRs) 6 66s ccc ct ee 
read(*,*)finl 
write(*,*) ‘Input Min ofeRRD imemneeGReors) .. os. sete s ee oe 
read(*,*)fin2 
write(*,*) 'Input MaxmoGmSR2 simmii .Ch1Oge. oc eee. + oe ee 
read(*,*) fin4 
write(*,*) ‘Input Max of web in IN . (Piotr... cee «ss. ee 
read(*, *) £inweb 
write(*,*)'Input increment for web in IN. (Rei O Sit crerctee so cs 
read(*,*)webstp 
write(4,20) fini, fin2, fin4, finweb, webstp 
else 

nmax=0 
write(*,*) 'Number of SpekeE?..«(E3).... ee ' 
read(*,*)nspoke 
ibb=-1 
write(*,*)'Do you want optional output? y-yes, n-no....... ' 
read(*,5)yes 
if(yes.eq.'Y'.or. yes. eq. 'y')then 

iprnt=-1 
else ae. 

iprnt=0 
endif 
write(4,10)nmax,nspoke, ibb, iprnt 
write(*,70) ; 
format(1lx,'Input initial port area, SQ IN..(F10.5)......... ') 
read(*,*)aport 
write(*,80) 
format(1x, 'Input initial web thickness, IN.(F10.5)......... ') 
read(*,*)web 
write(*,90) 
format(1x,'Input the perimeter factor..(F10.5).....eeeeeeee ") 
read(*,*)pfac 
kspoke=0.0 
write(*,60)rrs 
write(4,20)aport,web, pfac,kspoke,rrs 
write (*,*) "Input Ginitia bRRIMEBNCHES:. (Fle@n 5)... eee... o 
read(*,*)ril 
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write(*,*)'Input initial value of SRl (Default=0.25 IN.).... ' 
read(*,*)sril 
write(*,*)'Input initial value of SR2 (Default=0.25 IN.)..... ' 
read(*,*)sri2 
write(4, 20) cidyisril, sri2 
dell=0.0 
del2=0.0 
de13=0.0 
del4=0.0 
write(4,20)dell,del2,del13,del14 
cirtol=0.0 
pttol=0.0 | - 
pertol=0.0 a 
write (4,20)cirtol,pttol,pertol 
finl=0.0 
fin2=0.0 
fin4=0.0 
finweb=0.0 
webstp=0.0 
write(4,20)finl,fin2,fin4,finweb,webstp 
Wetec, *)einputeoqrain length, INCHES (F10.5) .cccccccesee ! 
read(*,*)lgrain 
write(*,*)'Input initial PHI2, DEGREES (F10. S)seseeseeeee, ! 
read(*,*)dphi2 
write(*,*)'Input initial PHI4, DEGREES (F10.5)...cccccceee ' 
read (*,*)dphi4 
write(4,20)lgrain,dphi2,dphni4 
write(*,*) 'Input # of steps in decimal format! 
write(*,*)'Input # of burnback steps for SR1>0......cecceee-' 
read(*,*)stepl z 
write(*,*)'Input # of burnback steps for PHI2>0 & SR1l=0.... ' 
read(*,*)step2 
write(*,*)'Input # of burnback steps for LSPOKE=0 & PHI2>0. ' 
read(*,*)step3 
write (*,*) 'Input # of burnback steps for ' 
Dement eG Etangular SPOKE OUPMNDACK. «ccc cc ecceecccceseee 
read(*,*)step4 
write(*, aoe tt of burnback steps for ' 
write(*,*)'remaining fuel DUINOUL. oer eeeeecesccececcreccecs : 
read(*,*)step5 
: write(4,20)stepl,step2,step3,step4,step5 

endif ey 
format (al) 
format (413) 
format (5f£10.5) + 
close(4, status='keep') ¥ 
return . ; : 
end 
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